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ABSTRACT 

A MICROPROCESSOR BASED ANTI-ALIASING FILTER 
FOR A PCM SYSTEM 

DANIEL CREIGHTON MORROW 
March 1983 

This project report describes the design and 
evaluation cf a microprocessor based digital filter- The 
filter was designed to investigate the feasibility of a 
digital replacement for the analog pre-sampling filters 
used in telemetry systems at the NASA Ames-Dryden Flight 
Research Facility (DFRF). The digital filter will utilize 
an Intel 2920 Analog Signal Processor (ASP) chip. Testing 
includes measurements of (a) the filter frequency response 
and (b) the filter signal resolution. The evaluation of 
the digital filter was made on the basis of circuit size, 
projected environmental stability and filter resolution. 
The 2920 based digital filter was found to meet or exceed 
the pre-sampling filter specifications for limited signal 
resolution applications. 
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CHAPTER 1 


INTRODUCTION 


A substantial portion of the research conducted at 
NASA's Ames-Dryden Flight Research Facility (DFRF) is 
flight oriented. Airborne test beds allow experiments to 
be conducted under conditions which cannot be approximated 
in wind tunnels or pressure chambers. Storing the data 
collected during an airborne experiment is a definite 
challenge in the space restricted flight environment. 
Often, the volume of data collected during flight tests 
make on-board storage of data, at best, impractical. 
Therefore, a pulse code modulation (PCM) 1 telemetry system 
is used to collect, multiplex and transmit data to the 
ground station for storage and analysis (see Figure 1-1). 
The ability of the PCM system to accurately measure and 
transmit data to the ground station is critical to the 
success of airborne experiments. 


^Ferrel G. Stremler, Introduction to Communication 
Systems , ed. David Cheng, Leonard A. Gould and Fred 
Manasse, (Ph i 1 1 i pi ne s : Add i son-Wesley Publishing Company, 
Inc., 1977) , pp. 355-366. 


1 
















The flight environment is full of electrical 
noise. The major contributors are the aircraft 400 Hz 
power supply, the avionics and the power plant. The pre- 
sampling filters (PSF) perform a twofold service for the 
system. In the process of attenuating noise signals in 
frequency bands above that of the data, these low pass 
filters act as anti-aliasing filters 2 ' 3,4 for the PCM 
system signal multiplexer. Without these filters, 
aliasing of high frequency noise into the data bands could 
seriously impair system accuracy. 

The Analog Pre-Sampling Filter 

The circuit diagram of the analog PSF is shown in 
Figure 1-2. If the resistors and capacitors of the 
circuit are matched, the transfer function of the circuit 
becomes: 

W 3 

H( S ) . r - 

s° + 2W C S" + 2W ( rs + 


Andreas Antoniou, Digital Filters: Analysis and 
Desian, (United States: McGraw-Hill , Inc., 1979) , pp. 13T- 
154.' 

^Lawrence R. Rabiner, Theory and Application of 
Digital Signal Processing , (New Jersey : Prentice Ha 1 1 , 
Inc., 1975), pp. 26-28. 

4 2920 Analog S ignal Processor Design Handbook , (Santa 
Clara: Intel Corporation, 1980), pp. 2-2, 2-3. 
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i 



R1 



n(s) = 

1 

As ?>s + Cs ■+ i 

where : 

A - R 1 *R 2 *R3 *C 1 *C2 *C3 
3 » RI* ( R2+R3 ) *C 1 *C2 
C = R1*(C1+C3) 


I 


Figure 1-2 Analog Pre- sampling Filter Circuit 
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where : 


w - 

1 

c 

RC 

R = R1 

= R2 = R3 

C = Cl 

= C2 = C3 


Thus, if the circuit is assembled using matched 
components, it provides a very efficient implementation of 
a third order Butterworth filter. The circuit requires 
few parts and occupies approximately two square inches on 
a circuit board. While the simplicity of this circuit 
allows a space efficient design, it also leads to 

sensitivity to environmental stresses. 

\ , 

\ . 

\ 

As the components age and go through the thermal 
extremes of flight, the values of the various circuit 
components tend to drift. The program TRCLC (see appendix 
E) was used to evaluate the magnitude and phase response 
of the circuit as the value of R1 varies. Figure 1-3 
shows the magnitude response curves when R1 is allowed to 
vary 10%. At the filter corner frequency (Wa = l Hz), the 
gain error associated with a 10% variation in R1 is 
1.2 dB. The distortions in the PSF transfer function lead 
to system errors in measured data. Fortunately, if the 
distorted transfer function is known, the impact on the 
final measurements can be minimized by post flight data 

i 
l 

1 

l 
i 


reduction . 
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The Digital Alternative 

A digital implementation of the PSF filters coaid 
provide a significant improvement to the existing PCM 
system. The transfer functions of digital filters are 
primarily dependent on the program which drives the 
processor, rather than the values of the components which 
form the circuit. Therefore, the digital filter is 
expected to exhibit improved transfer function stability. 
By the same argument, improved filter repeatability is 
also expected. The potential also exists for reducing the 
physical size of the filter. 

This project report describes the design and 
evaluation of a real-time microprocessor based digital 
PSF. The goal of the project was to investigate the 
properties of a digital filter which approximates the 
analog PSF circuit. The emphasis of the design is on size 
efficiency and transfer function stability, not on the 
implementation of an ideal Butterworth filter. The 
processor utilized to implement the digital filter is the 
Intel 2920 Analog Signal Processor (ASP). This processor 
was selected because of its unique design and the 
availability of support tools. 



CHAPTER 2 


DIGITAL FILTER DESIGN 


Filter Design Specifications 

The pre-sampling filter is an important part of 
the PCM system. The filter attenuates high frequency 
noise in sensor signals before the PCM system samples the 
data for transmission. The high frequency noise must be 
removed to avoid aliasing. The performance of the 
existing analog filters is degraded by their sensitivity 
to environmental stresses. The digital PSF is to 
approximate the response characteristics of the analog 
circuit. The nominal specifications of this filter are: 


1) The filter transfer function: 


-G 

H ( s) 

S , S St 

( ]3 + 3 [ ]2 + 3 [ ] t 1 

w c W c W c 

where: s = jW 

W c = Corner frequency 

G = Ga i n 

2) A balanced differential input. 

3) Input impedance: 100k ohms minumum. 


8 
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4) Output impedance: 20k ohms maximum. 

5) Nine bit {eight bit plus sign) signal 
resolution . 

6) Four signal inputs per 2920 chip. 

The digital filter design process will be treated 
in four major steps: 

1) Filter hardware design 

2) Choosing a sample frequency 

3) Development of the Z domain transfer function 

4) Implementation of the transfer function 

The following is a discussion of the engineering 
process for designing a digital replacement for the analog 
third order Butterworth PSF. 

The 2920 Support Hardware Design 

The 2920 ASP requires several types of external 
hardware support to meet the PSF nominal specifications. 
The 2920 input channels require buffering to meet 
impedance requirements and to interface with the 
differential input signals. The output signals require 
filtering to remove noise from the D/A converters. 

Figure 2-1 shows the functional blocks necessary 
to interface a 2920 based filter with the existing PCM 


Ili f ferent i.al Input s 

VV VV VV YV 


10 
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Figure 


-1 Digital P3F Functional 3 locks 


f'SF. Outputs 
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system. The Signal Conditioning Buffers (SCB) are 

/ 

required to meet the nominal filter input specifications. 
The SCB circuits raise the input impedance of the filter, 
convert the differential inputs to single line, scale the 
signal to magnitudes acceptable for the 2920 ADC's and act 
as anti-aliasing filters for the 2920 by attenuating high 
frequency signals. The 2920 output signals are passed 
through reconstruction filters to attenuate noise from the 
DAC outputs. The reference voltage will be generated 
local to the processor to reduce reference voltage noise. 
The circuit shown in Figure B-l is one possible 
implementation of these functional blocks. Further 
information on 2920 support hardware design is contained 
in Appendix B. 

Choosing Sample Frequency 

Choosing the digital filter sampling frequency 
(F s ) is a fundamental part of the digital filter software 
design. As shown in Figure 2-2, a broad range of sample 
frequencies is available for 2920 filter design. The 
acceptable values of F s are band limited by the Nyquist 
frequency-^, noise bands and the speed of the digital 
processor. The minimum acceptable sample rate is 


^■Ferrel G. Stremler, Introduction to Com munication 
Systems, ed. David Cheng, Leonard A. G o u l"3 and Fred 
Manasse, (Phill ipines: Add ison-Wesley Publishing Company, 
Inc., 1977) , pp. 112-116. 
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ORIGINAL 
OF POOR 


pkQB'- & 




Indicates regions of unacceptable sanding 
freguenc i es 


Fc = 10 Hz 

Fsmax =34.7 KHz for 2920-16 
= 26.0 r'.Hz for 2920-13 



Figure 2-2 Acceptable Sampling Rates for a 2520 Based 
' 1 0 Hz Filter 

• ( . 



determined by the Nyquist criteria applied to the highest 
signal frequency of interest. This frequency corresponds 
to the point at which the transfer function of the filter 
is down 54 dB (the maximum resolution of the 2920 inputs 
and outputs). This occurs approximately a decade above 
the corner frequency of the third order Butterworth 
transfer function (10 F c ). The Nyquist criteria^ states 
that the minimum acceptable sample frequency (f’gm'in) .-M.- 
twice the -^54 dB frequency or: 


F cm , = 20 F _ 

sntin c 


where: F = corner frequency of the f i Iter 


F s should be chosen such that it will not cause 
aliasing of the 400 Hz power supply noise into the 
passband of the digital filter. The frequency bands which 
are subject to aliasing are. centered around multiples of 
F , and are bounded by: 


N F s - 10 F c 
and 

N F s + 10 F c 

where N = 1, 2, 3, 4, . . . 


^Streraler, pp. 112-116. 
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Any noise occurring in these bands is frequency 
translated by the sampling process into the passband. 
Therefore, to avoid frequency aliasing of the 400 Hz power 
supply noise signal into the passband, F g must be chosen 
such that: 


N F c - 10 F > 400 

o C 

or 

N F s + 10 F c < 400 

Sample rates which do not meet this condition are 
considered unacceptable for this application. The maximum 
sample rate ( F smax ) is established by the maximum clock 
frequency of the 2920 processor. This rate as a function 
of crystal frequency is: 


smax 


cry 


N 


ps 


where: F = crystal frequency of the processor 
cry 

F cry < 6 * 6E6 for 2920-16 processor 

F < 5.0E6 for 2920-18 processor 
cry 

Np S = number of program steps in the filter code 


The program cycle rate (Fp Cr ) is the frequency at 
which the program space of the 2920 is cycled. The value 
of F s is also constrained to be a factor or multiple of 

F pcr* F s equals F pcr 


if a signal sample is processed 
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every program pass. F g may be made less than Fp Cr by a 
factor N, by enabling the filter input sequence on every 
Nth program pass. This technique is discussed in greater 
detail in the 2920 Ana log Signal Processor Design 
Handbook . F s may be up to four times greater than Fp cr by 
linking duplicate filter code blocks in the 192 step 
program memory. 


In Summary, the acceptable values of Fs are 
limited by the following constraints: 


1 ) 


2 ) 

3) 


F = F * 
s cpr 


M 


1 

where: (M=4,3,2,l, — 

2 


F_ > 20 F_ 
s c 

N F s - 10 F c > 400 
or 



1 


4 


) 


N F s + 10 F c < 400 

where': Fp cr = program cycle rate 

F 

cry 
4 N 

ps 

and: F = processor crystal frequency, 

c t y 

F c = filter corner frequency, 

Np S = number of program steps in the 


filter code. 
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Development o£ the Z Domain Transfer Function 

The development of filter software which will 

approximate the characteristics of the analog PSF requires 

that a Z domain equivalent of H(s) be developed. The type 

of transformation process which will yield the best 

results depends on the form of H(s) and the target filter 

structure. For this application, which requires a space 

efficient implementation, Infinite Impulse Response (HR) 

transform techniques will yield the best results. IIR 

filters approximate analog circuit transfer functions more 

efficiently than Finite Impulse Response (FIR) filters. 

Further, studies show that the Bilinear Transform is the 

preferred IIR transform technique for Butterworth transfer 
■3 

functions . 

The B i linear Transform 

The Bilinear Transform^ technique is an 
approximation method which converts an analog transfer 
function H(s) to a d iscrete-time transfer function H(z). 
This method allows any valid analog transfer function H(s) 


Bernard Gold and Lawrence R. Rabiner, Theory and 
Appl i cat ion o f Digital Signal Processing , (New Jersey: 
Prentice-Hall, Inc., 1975), pp. 252-257. 

‘‘Gold and Rabiner, pp. 219-224. 



to be expressed as a function of the discrete operator z 
using the equation: 
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H ( z ) = H (s) 


s 


2 F, 


z - 1 

C ] 

z + 1 


The Bilinear Transform is an algebraic 
substitution of a function of the discrete variable z for 
the laplacean variable s . The derivation of this 

technique is based on a trapezoidal approximation of an 
integral, thus the Bilinear Transform technique is not an 
exact method . The frequency response of H(z) is equal to 
H(s) if and only if : 

W d 

W a = 2 F Tan [ ] 

2 F 
* s 

where: W_ = frequency variable (analog) 

Q 

= frequency variable (discrete) 


This equation implies that the frequency response 
of the Bilinear Transform based discrete time transfer 
function will be frequency shifted with respect to H(s) as 


^Andreas Antoniou, Digital Filters: Analysis and 
Design, (United States: McGraw-Hill, Inc., 1979), pp.178- 
179. 
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a function of input frequency. This effect is known as 
" w.arping" 6 . 

The warping error is negligible for input 
frequencies much less than F s . The warping error 
increases as the frequency variable approaches one-half 
of the sample rate. The error may be minimized in the 
pass band by applying the pre-warping factor (L w ) during 
the transform process. L w is defined as: 


2F S Tan [ ] 


where: W c = filter corner frequency. 


L shifts the frequency response of the digital 
transfer function so that the amplitude response of the 
digital transform is exactly equal to the amplitude 
response of the analog transfer function at W c . Applying 
the Bilinear Transform to the transfer function of the 
general third order Butterworth pre-sampling filter shown 


^Antoniou, pp. 182-185 



in Figure 1-2 we find: 


H ( s ) = 


S'? s ? si 

[ — r+ 3 [ — ] 2 + 3 [ — ] x + i 


w. 


W, 


W, 


s=2F s L w [ — 1 


.z-1. 

z+1 


yields the z transform; 


H ( Z ) = 


where: K 
K 
K 
K 



lz 

3 + 

3z 2 

+ 

3z 

+ 

1 



K o 

z 3 

+ 

Z 2 

+ K 2 

z 

+ 

K 3 



: o = 

8 

R ^ 
R sc 

+ 8 

r 2 

sc 


+ 

4 R sc 

+ 

1 

1 

24 

R 3 
sc 

8 

r 2 

sc 


+ 

4 R sc 

+ 

3 

2 = 

24 

R ^ 
R sc 

8 

R2 

sc 


- 

4 R sc 

■f 

3 

3 = 

-8 

R ^ 
H sc 

+ 8 

R 2 
R sc 


- 

4 R sc 

+ 

X 


R s „ = sample rate to corner frequency ratio 


r * F 
L w r s 


W c 

2 Tan [ ] 

2F„ 


z = discrete-time operator 


The analog to digital transfer function 
transformation for the Butterworth filter is complete once 
the values of K Q through K-^ have been calculated. The 
results of the transformation are compared to the analog 
response in Figure 2-3. The frequency response of the 
H ( z) was evaluated by the program MAGCK (see appendix E). 







Normalized Frequency in llz (I.og Scale) 
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The magnitude response of the transforms is normalized to 
0.8 volts/volt to avoid the non-linearities in the 2920- 
16D ADC's and DAC's. This normalization convention is 
carried throughout this investigation. 


Implementation of H ( z) 

Implementing H(z) with a 2920 ASP requires that 
the discrete transform be converted into a sequence of 
instructions that the 2920 can execute. The structure of 
the filter algorithms is an important consideration given 
the program length restrictions of the 2920. Realizing 
four filters with one 2920 implies that the length of each 
individual filter algorithm is limited to forty-eight 
program steps. 

Benchmark programs were written based on the 
ladder, parallel, direct and cascade structures to 
evaluate the relative efficiencies of the 2920 algorithms. 
The benchmark filters implement 300 Hz PSF's with a 
minimum signal resolution of nine bits. The parallel and 
ladder structures did not perform well in this 
application. Both these structures failed to meet the 
signal resolution requirements. The direct and cascade 
forms were found to be better suited to the 2920 


architecture . 
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Figure 2-4 shows the flow diagrams of the cascade 
and direct form-2 algorithm structures 7 . The Z“^ term 
represents a memory storage element where data is time 
delayed one sample. The sections with discrete feedback 
are called recursive sections, the sections without 
feedback are called non-recursive sections. The recursive 
sections implement the denominator (poles) of H(z), and 
the non-recursive sections implement the numerator 
(zeros) . 


The direct form-2 structure implements H(z) with 
fewer multiplications and memory variables than the 
cascade form, but the resolution requirements impact the 
length of the multiplication algorithms to a greater 
degree than in the cascade form. The direct form requires 
three full accuracy recursive multiplications and four 
non- r ec ur s i ve multiplications. The cascade form also 
requires four full accuracy recursive multiplications, but 
the gain of the recursive section is divided into two 
parts, thus reducing the accuracy requirements of these 
multiplications. Further, the gain of the two non- 

recursive sections can be distributed so that three of the 
multiplications are one instruction long. This 
effectively reduces the number of non-recursive 


7 Antoniou, p.76. 
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multiplications in the cascade structure to two. 

The third order recursive section of the direct 
form is subject to oscillations caused by power up 
conditions. Under normal circumstances stability is not a 
problem for these digital Butterworth filters. 
Occasionally power up conditions leave random values in 
the delay memories which cause ALU overflows during 
multiplications. Recursive sections feed the overflow 
resultants back into the algorithm. Under certain 
conditions the overflow resultants cause new overflow 
errors. When this occurs the recursive section goes into 
an oscillatory mode. The instability is observed only in 
odd orders of recursive sections. This unstable mode can 
be prevented by software which detects the mode and damps 
the oscillations. The additional instructions of the 
software fix further reduce the efficiency of the direct 
structure algorithms. 

A comparison of the relative efficiencies of the 
benchmark programs shows that the cascade structure is 
best suited for the PSF application. 


The Cascade Difference Equation 

The cascade form of H(z) may be derived directly 
from the results of the Bilinear Transformation. Applying 
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the Bilinear Transform to H(s) it was found that: 


Y(z) 

lz 3 

+ 3z 2 + 3z 

+ 

1 

X(z) 

Kq z 3 

+ z 2 + K 2 z 

+ 

K 3 


The cascade form of the discrete transfer function 
is derived by expressing H(z) as the product of second 
order transfer functions. In this application, H(z) will 
be expressed as the product of two transfer functions, 
G^Cz) and G 2 (z). The transfer functions G-^z) and G 2 (z) 
are defined as: 

Y(z) 

H(z) = = G (z) * G 2 (Z) 

X(z) 


where : 


G 1 (z) = 


ana : 


G 0 (z) = 


G' (z) 

1+ 2 z 1 + 1 z 2 

X(z) 

-1 -7 

l-d 1 z J ’-d 2 z 

Y(z) 

a o + a i Z_1 

G' (z) 

1 - d 3 z 1 

^ are defined such that: 

1. 

N(z) = — 

- [ 1 + 3 z -1 + 3 z -2 + 1 z 

K 0 


= [ 1 

+ 2 z _1 + 1 z -2 ] * [ a Q + 


.-1 
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3 r.c d^, d 2 » d^ are defined such that: 

D(z)= 1 - - b 2 z -2 - b 3 z ~ ^ 

= [ 1 - d 2 z' 1 - d 2 z" 2 ] * [ 1 -d 3 z" 1 ] 

with: b^ = -K^/Kq 

b 2 = ~ K 2 / ' K 0 

b 3 = -k 3 /k 0 

Gain Scaling in The Cascade Algor i thm 

If the flow diagram of Figure 2-4a was programmed 
as shown into a 2920, the output Y(nT) would be distorted 
by ALU overflow errors. The two’s complement fixed point 
2920 ALU overflows if the result of any calculation 
exceeds one (some overflow protection is provided, see the 
2920 Assembly Language Manual for details). To avoid ALU 
overflow errors the signal should be scaled as it is 
processed through the filter algorithm. 

The magnitude of a section gain determines the 
amount of signal scaling required to prevent overflow 
errors. A gain scale factor, C i , i s carefully introduced 
into each section such that the gain of the entire filter 
does not change. This requires that: 

i = C 1 * C 2 * c 3 * 

where: i = section number 



OR56SKSM. !S 
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A further restriction on the value of is made 
- 0 insure that the scale factors do not significantly 
lengthen the instruction code. is restricted to be a 

-ower of two. This allows to be implemented with one 
instruction. The binary scale factor ( S F ^ ) is the 
exponent of two which determines C^. The values of SF^ 
ere limited to integers between negative thirteen and 
positive two by the instruction set of the 2920. The two 
restrictions on the allowable values of C,- are: 


1) 


2 ) 



where SF i = -13,-12,-11,-10, 0,1,2 

! = C 1 * C 2 * c 3 * Cj 


Signal scaling prevents ALU overflows, but it also 
increases the accuracy requirements in the algorithms, 
which increases the required accuracy of the 
multiplications. The optimum magnitude of will be the 
largest valid power of two which prevents ALU overflows 
during signal processing. This implies that: 

. C i A i (l+k) < 1 

: 1 

i * Solving for C.- : C - < 

i ! 1 1 Aj { l+k ) 

i 

! i • 

i i 

i j 

I J where: Cj = gain scale factor 



i 
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and: A^ = gain of the section 

I 

= G(z) I 

I z=l 

with G { z ) = the section transfer function 

k = percent overshoot of the section 
transfer function 

Solving this equation for the binary scale factor SF^ : 

Ln + Ln (1+k) 

Ln 2 

The binary scale factor may be determined from 
this equation by finding the largest integer value of SF^ 
which will satisfy the equation. 

The method in which gain scaling is implemented is 
different for the recursive and non-recursive sections. 
In recursive sections, the scale factor is implemented as 
a binary shift applied to the section input. The 
coefficients are not modified. In non-recursive sections 
the coefficients and the scale factor are multiplied 
together to form new non-recursive section coefficients. 

Once the gain scale factors are properly placed in 
the flow diagram (see Figure 2-5), it is ready to be 




X ( n T ) 



Figure 2-5 The Flow Diagram of the Digital Filter Test Algorithms 
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converted to 2920 instruction code. C i and C Q buffer the 
magnitude of the signal input and the algorithm output. 


Accuracy o f Multiplications 

The flow diagram in Figure 2-5 contains five full 
accuracy multiplications. 2920 program memory constraints 
limit the number of instructions in each filter algorithm 
to forty -eight steps. An average 24 bit multiplication 
requires eight program steps. Five 24 bit multiplications 
would consume forty of the forty-eight available program 
steps. Clearly, if all the filter functions are to be 
implemented within the forty-eight instruction step limit, 
the multiplication algorithms must occupy less space. The 
multiplication algorithms can be made more space efficient 
by decreasing the accuracy of the multiplications. The 
amount of multiplication accuracy required to preserve the 
signal resolution, as it is processed through a section, 
is a function of the resolution of the input signal and 
gain of the section. 

The binary resolution of a fixed point digital 
signal is defined as the number of bits in the signal 
which contain accurate data. The original accuracy of the 
sampled data is limited to eight bits plus sign by the 
resolution of the A/D converters. The binary resolution 
of a signal ( R ^ ) at the input to section i can be 



integer which 
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approximated by finding 
satisfies the equation: 


the largest 


R i < R s - SF i 


Ln A 


Ln 2 


where: R. = Binary resolution of input to section i 
( number of bits) , 


R s = Binary resolution of signal 
(number of bits) , 

A = Product of the filter sections prior 
to section i. 

To establish the relationship between the input 
resolution, section gain and multiplication accuracy, the 
effects of finite signal accuracy in the section 
algorithms must be established. 


The value of a finite accuracy signal may be 
expressed as the value of the infinite accuracy signal 
plus an error term (E rs ). 

X(nT) f = X(nT) in + E rs 

where: X(nT)^ = value of finite accuracy signal 

X(nT)in = value of infinite accuracy signal 
E rs = finite signal accuracy error 
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Recursive section difference equations have the form: 

X(nT) in = Y(nT) in + B x X(nT-T) in + B 2 X(nT-2T)- n ... 
...Bj X(nT-jT) in 

where: Y(nT) = discrete input to the section 
X(nT-T) = section output 


Then, assuming finite signal accuracy, the 
difference equation becomes: 

X(nT) f = *(nT) in + E rs + (B^X (nT-T) - n ) + (B^E^) 

+ (B 2 *X(nT-2T) in ) + (% 2 * E rs ) + ... 

... + (Bj*X(nT-jT) in ) + (B j *E rs ) 

or: X(nT) f = X(nT) in + E rs + (B 1 *E rs ) + ( B 2 * E rs ) + •** 

... + (B j *E rs ) 

The finite signal error E rs is bounded by the 
signal resolution: 


E rs ! < 


-Rj 


Therefore, the maximum error contributed to a 
.signal by the multiplications in a recursive section may 
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be expressed as: 


F 

rs i rsec 

= E _ + 

B 

i * E 

+ B_ + 

E 

B. * E „ 

r s 


1 rs 

2 

rs 

j r 


= E * 

rs 

( 

1 + B 1 

+ B 2 + 

B 2 ... ) 


< 2 


-Ri 


* 


1 


A i 


where: E rs | rsec = maximum section output 
error due to finite accuracy sig- 
nal input for recursive section. 


Similarly, for the non-recursive section it can be 
shown that: 


c* = n * f + r> * r + n * p 

"rs I nrsec u 0 "rs U 1 h rs u 2 “rs 


D - * E 

u j c rs 


* E rs * < D o + D 1 + D 2 + 


• D.) 



where: E rs | nrsec = maximum section output 
error due to finite accuracy sig- 
nal input for recursive section. 

and: A^ = gain of section i. 


The effects of finite accuracy multiplications are 
determined by a similar argument. The multiplication 
are expressed as the sum of the infinite 


constants 


accuracy constant (B-) plus an error (E ) 

3 r n 


B. = 3 • + E 


3f 


3 in r j 


where: B 


3 f 


B 3 in 


r 3 


the value of the finite accuracy 
constant, 

The value of the infinite accuracy 
constant , 

finite constant error for Bj. 


Again, recursive section difference equations have 
the form: 

X(nT) in = Y(nT) in + (B 1 .* X(nT-T) in ) + (B 2 .* X (nT-T ) in ) • • • 

... + (B . * X(nT-jT) ) 

3 in 


The results of the finite accuracy multiplications 
for recursive sections may then expressed as: 

X(nT) f = Y(nT) in + (B, + E_ ) * X(nT-T) + ... 

r in i ln r^ 

. . . + (B, + E_ . ) * X (nT-jT) 

3 i n r j 

= X(nT) f + E * X (nT-T ) + E * X(nT-2T) +... 

... + E * X (nT-jT) 
r 3 

X(nT) is bounded by the fixed point memory such that: 


I X(nT) I < 

I I 


M 
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Also , 

bounded by the 


the individual 
accuracy of the 


multiplication errors a 
multiplication algorithm: 


e 


I I 

I E . I '< 2 

I r 3 I 




where: R m = accuracy of multiplications (as- 
sumed to be uniform in section i) 


Therefore the the maximum output error associated 
with finite accuracy multiplications in recursive sections 
is : 


i < N * 2 

rm I r sec m 


-R 


m 


1 . 


A i 


where: E rm | rsec = maximum output error due to 

finite multiplications in a 
recursive sections. 


For the non-recursive sections, the error associ- 
ated with finite multiplications is: 


rm | nrsec 


< N. 


-R 


ra 


ra 


* A . 

A i 


where: E rm | nrsec = maximum output error due to 

finite multiplications in a 
non-recursive sections. 
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The finite accuracy multiplications will not 
degrade the section performance if the error associated 
with the multiplications is less than the maximum output 
signal resolution. Un-degraded section performance 
requires that the errors due to the finite multiplications 
are less than errors due to the finite signal accuracy. 
The inequalities: 


F < r 

rm | rsec rs | rsec 

and 

F < p 

rmlnrsec rsjnrsec 

state this condition. Solving these inequalities for the 
required binary multiplication accuracy yields the same 
relationship between multiplication accuracy, number of 
multiplications and input signal resolution for recursive 
and non-recursive sections. 

These equations reduce to the inequality: 

Ln N m 

R m > R i + I 3 

Ln 2 

where: = required binary accuracy of 

section multiplication 
= input signal resolution 
N m = number of multiplications in 


the section 
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The minimum required accuracy of the 
multiplications in a section may be determined by finding 
the smallest integer which will satisfy this equation. 

Converting the final flow diagram to 2920 instruc- 

\ 

tion code is a straight forward task. The coefficient 
multiplications must be coded into efficient algorithms 
using the section Rm as a standard for minimum accuracy. 
Figure 2-6 shows a breakdown of the digital and analog 
function placement in the forty-eight step program block. 
The code is arranged to minimize the phase delay 
associated with numerical processing. The code block 
shown was used in the development of the PSFSS software. 
The A/D conversion is not complete until late in the 
processing, so all the multiplications which do not 
require the input are performed first. The final portion 
of the code block performs the sample rate division 
process and manages the time delayed variable storage. 










CHAPTER 3 


RESULTS AND MEASUREMENTS 

The evaluation of the 2920 based PSF includes 
several areas of its design. The filter size, accuracy 
and limitations are all important considerations. The 
issue of filter size is investigated in Appendix B. The 
test results presented in this section demonstrate filter 
accuracy and provide an indication of the limitations of 
the 2920 ASP. 


The final space requirements of a 2920 based PSF 
are largely decided by the complexity of the 2920 support 
hardware. The 2920 occupies 1.5 square - inches of the 
estimated 10 square inches occupied by the PSF circuit 
shown in Figure B-i. The remainder of the space is 
consumed by the various 2920 support and PCM interface 
circuits. The circuit presented in Appendix B was 
designed using available components in order to provide a 
rough estimate of the size requirements of a 2920 PSF 
which meets all the input and output specifications. The 
design is not space optimal. Several special purpose 
linear chips available could be used to reduce the size of 
the circuit considerably. 
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The accuracy of the 2920 processor and PSF 
algorithms was determined by measuring filter transfer 
functions from a 2920 PSF test circuit. Three basic PSF 
algorithms were written for the 2920, and from these, 
using sample frequency division techniques, fourteen 
filter algorithms were developed. These algorithms form 
the P3FSS library presented in Appendix D. 


Test C i rcu i t ry and Transfer Function Measurements 

A block diagram of the test set which measured the 
transfer functions is shown in Figure 3-1. The filter 
input signal is provided by a noise source from the HP 
5 4 4 1 0 A A/D Converter. The output of the filter is 
processed by the HP 54410A A/D Converter and a HP 54470B 
Digital Filter. The digitized output signal is then 
processed by the HP 5420A Digital Signal Analyzer to 
determine the filter transfer function. A circuit diagram 
of the 29 20 and test support hardware is shown in Figure 
3 - 2 . The input SCB's and local reference voltage 
generator were not included in the test set to eliminate 
the impact of these circuits on the measurements. The 
power supply voltages, reference voltage and bread board 
for the test circuit were supplied by an EL Instruments 
ELITE 2 Circuit Design Test System. The reconstruction 
filters on the 2920 outputs were required to minimize 
interaction between the 2920 DAC's and the sampler in the 
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Figure 3-1 Filter Test Set 
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HP A/D converter. Without these filters the high 
frequency noise added by the D/A conversion could cause 
aliasing errors in the measurements. The transfer 
function of the test circuit was measured using a 2920 
test program which passes the sampled input directly to 
the output without modification (other than a D-C offset 
error compensation). The test circuit transfer function 
(Figure 3-3) illustrates several effects caused by the 
sampling process^. The magnitude response of the circuit 
shows the magnitude of the output as a function of the 
input frequency. The frequency of the output signal is 
always less than the sample rate, regardless of the input 
frequency. The magnitude minimum and 180 degree phase 
shift observed at multiples of the sample rate are 
characteristic of the test set. The slight roll off of 
the magnitude response is caused by signal skewing during 
the sample capacitor charging cycle. 

2920 Version D Software Fixes 

The D version 2920 has several glitches which 

o 

affect the processor operation . Intel has released new 


•*■2920 Analog Signal Processor Design H andbook , (Santa 
Clararlntel Corporation, 1980), pp. 2-4. 


2 Carmel Vernia, "Signal Processing Products Status", 
Signal Processing News , April 1982, p. 2. 






46 


versions of the 2920 which are free from these operational 
problems'*. Most of the operational 'glitches' had no 
impact on the PSF filter operation; others require 
software fixes to insure proper operation. These fixes 
appear in the input sequence of all the filter algorithms 
written for this investigation. The SUB DAR,KM2,R00 ,CND6 
and LDA T,T,R00,CND4 instructions are all version D fixes. 
These instructions occupy four program steps which could 
otherwise be used for multiplications. Separate fixes are 
used to compensate for signal conversion errors. 

A test algorithm which samples the analog input, 
converts it to a digital value, then back to analog was 
used to measure the 2 9 20 conversion process 
characteristics. Figure 3-4 shows the results of the 
measurement. The linear regions for the output converters 
are from -.95 to 0 volts, and 0 to .9 volts. Beyond these 
regions, the converters saturate. The gain of the PSF 
filters has been normalized to .8 volts/volt to prevent 
DAC saturation. The conversion processes also introduce a 
negative DC offset. A typical offset can be seen in the 
linearity curve shown in Figure 3-4. Intel states this 
offset is typically -16 mV, -70 mV maximum. To partially 
correct for the offset, a DC value of .062 is added to the 


^Carmel Vernia, ”2920 Status", Signal P rocess i ng 
News, January 1982, pp. 6-12. 



Vout ( vol ts ) 


47 



Figure 3-4 2920 Convertor Linearity 
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filter output before the D/A conversion. These PSF fixes 
appear in all the filter algorithms presented in this 
investigation. Five digital and four analog instructions 
in each forty-eight instruction block are dedicated to the 
version D fixes. 

Test Algor i thms 

The test algorithms presented represent the 
extremes of the filter applications. The test set 

includes 300 r 80, 70 and 10 Hz filters. All the 
algorithms we re written using the cascade structure of 
Figure 2-5, and are part of the PSFSS library (see 

Appendix D) . The 300, 80 and 70 Hz algorithms sample at 

4. IK Hz. The 10 Hz algorithm is actually the 80 Hz 

algorithm using the sample rate division techniques to 
reduce the sample rate and corner frequency by a factor of 
eight. The 300 Hz and 10 Hz filters represent the 

extremes of the PSF corner frequencies used by NASA DFRF. 
The 70 Hz filter represents the approximate lower limit of 
the filter corner frequencies without using sample rate 
division techniques. 

Magnitude and Phase Response Curves 

The transfer functions of the algorithms as 
measured and plotted by the test set are shown in Figures 
3-5, 3-6, 3-7 and 3-8. The expected transfer functions 
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are shown in Figures 3-9, 3-10, 3-11, 3-12. The measured 
transfer functions correlate closely with the predicted 
transfer functions. The 2920 filter resolution is 
demonstrated in Figure 3-13. This plot shows the high 
frequency response of the 10 Hz filter. The filter rolls 
off till the magnitude of the output signal is less than 
the resolution of the DAC's. The result is the flat area 
in the stop band above 100 Hz. This plot also shows the 
effects of sample aliasing. A slight rise in the 
magnitude response may be observed at 200 Hz. This effect 
is shown more clearly in Figure 3-14. In this plot, the 
frequency scale of the 10 Hz filter transfer function has 
been expanded to show the pass bands centered at multiples 
of the sample rate. The peaks are caused by ADC frequency 
aliasing of the analog input signals. 

Filter Repeatabi 1 i ty 

A demonstration of filter repeatability is shown 
in Figures 3-15 and 3-16. These plots are the transfer 
functions of two 70 Hz filters implemented on different 
2920 chips. The frequency responses of these filters are 
for all practical purposes identical. 
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Figure 3-9 Calculated Magnitude Response of the 300 Kz 
Digital Filter 
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Calculated Magnitude Response of the 70 K 
Filter 
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Calculated Magnitude Response of the 10 Hz 
Filter . 
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Figure 3-13 Magnitude Response of the 10 Hz Filter 
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CHAPTER 4 


CONCLUSION 


Summary 

The goal of this investigation was to evaluate the 
feasibility of replacing analog PCM PSF's with a digital 
equivalent based on Intel's 2920-D-18 Analog Signal 
Processor. The focus of the investigation was on the 
operational characteristics and limitations of the 2920 
processor functioning as a PSF. 

The design procedure established for this 
investigation starts with the analog PSF transfer function 
H(s). The Bilinear Transform Technique is applied to K(s) 
to derive the digital transfer function H(z). H(z) is 
then factored into its cascade form. Scale factors are 
assigned on the basis of individual filter section gain to 
prevent ALU overflows during processing. The resulting 
modified cascade transfer functions are coded into 2920 
assembly language. These filter algorithms are then 
compiled and loaded into the 2920 EPROM memory. The 
filter algorithms were ground tested in a 2920 test 
circuit to determine the operating characteristics of the 
filter . 
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The test results presented in the previous chapter 
demonstrate that the 2920 is capable of performing the PSF 
function for limited resolution signals. The maximum 
signal resolution which the 2920 can process without some 
loss of signal accuracy is determined by the width of the 
2920 A/D and D/A converters. This constraint limits the 
resolution of 2920 output signals to nine bits. There are 
several design 'tricks' which extend the resolution of the. 
2920, but the techniques are too hardware and software 
intensive to be used in this application. 

The transfer functions of the digital filters roll 
off more sharply than the analog filters, and have 
excellent repeatability characteristics as the filter 
algorithm is implemented in different chips. The 
increased roll off rate is due to the "warping effect" of 
the Bilinear Transform. The frequency scale distortion in 
the H(s) to H(z) mapping is desirable in this application. 
The increased noise rejection in the transition band is 
realized without disturbing the smoothness of the 
Sutterworth shape. The transportability of filter 
algorithms is a major advantage of the digital filter. 

The filter algorithms were written in such a form 
that a computer program can handle all the code linking 
and compiling operations normally associated with digital 
programs. The filter transfer functions are independent 
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of which 2920 chip or where in memory it is implemented. 
Therefore, the algorithm development process can be 
performed in a uniform manner for any 2920 chip. The 
PSFSS programs are a result of this characteristic of 
digital processing. The PSFSS tool simplifies the 
installation of a 2920 PSF modification to a five minute 
operation. The only information needed to specify a PSF 
is the desired corner frequency. The power in this type 
of tool lies both in the simplicity and the flexibility. 
The PSFSS software manages signal processing modules. The 
function of the modules is not limited to PSF filters. 
Modules could be included which perform any desired 
function (waveform generation, linearization, etc.). This 
flexibility could allow many of the tasks normally 
associated with ground processing to be performed real 
time in flight. 

The size of the 2920 based PSF circuitry is 
estimated to be 2.25 square inches per filter. This is 
comparable to the size of the analog circuit. This 
estimate is based on the circuit shown in Figure B-l, and 
the implementation of four filters per 2920 PSF circuit. 
The potential exists for significantly reducing the size 


of the circuit. 



2920 Version J 


The major design limitation encountered in this 

investigation was the length of the 2920 program space. 

Implementing four filters per chip limited the number of 

instructions for each filter algorithm to forty-eight. 

This design constraint becomes more and more imposing as 

the ratio R sc (f’ s /W c ) increases. As R sc increases, the 

accuracy requirements of the multiplication .algorithms 

also increase. More accurate multiplications require more 

instruction steps. If R g „ is large, there simply is not 

enough room to implement the algorithm within the forty- 

eight program step limit. The 70 Hz filter in the PSFSS 

library has a Rsc equal to 9.5 Hz/Rad. The upper limit of 

R which can be imolemented within the size and 
s c 

resolution constraints is approximately ten for the 
2920-D. The filter algorithms with corner frequencies 
below 70 Hz (R sc > 10.0) were created using the sample 
rate division techniques. The value of R sc for these 
filters was too large to allow coding a filter which used 
the base sample rate (4166.7 Hz). Intel has released a 
newer version of the 2920 which does not require the 
software fixes required by the version D 2920 used in this 
investigation. The I/O sequences version J processor have 
been fixed thus increasing the number of instructions 
available for calculations. The impact of the reduced 
algorithm overhead of the version J is to increase the 
flexibility of algorithm design. 
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Areas for Further Study 

This investigation demonstrated the 
characteristics of the 2920 PSF transfer functions. But 
the effect of environmental stresses on the critical 2920 
support circuits has not been established. Further study 
and tests are required to determine an optimum design for 
the support hardware and the impact of environmental 
stress on the 2920 PSF circuitry. 

The function of the analog PSF in the PCM system 
is to remove high frequency noise in signals before 
sampling. The introducing of a digital equivalent for the 
analog filters creates an interesting anomaly in the 
system. Analog signals are processed by the 2920 circuit 
SCB’s to remove any high frequency noise, then sampled by 
the 2920. The 2920 processes the digitized signal, then 
converts the signal back to analog so it can be sampled by 
the PCM system. The PCM sampler has a digital PSF (the 
2920) which has its own PSF (the SCB's). Further, the 
signal is converted from an analog to a digital twice. 
The extra PSF and conversion could be eliminated by using 
the SC B circuits as PSF's for the PCM sampler and 
performing the digital filter processing function on the 
output of the PCM sampler. Further study is needed to 
determine the feasibility of eliminating the extra PSF 
circuit and A/D, D/A conversion process present in the 
2920 based PSF scheme. 
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APPENDIX A 


TRANSFER FUNCTION GRAPHS 


Appendix A contains the measured transfer 
functions of the 2920 algorithms contained in the PSFSS 
filter library- The transfer functions were measured and 
plotted using the test set described in Chapter 3. 
Figures A-l through A-14 are the magnitude and phase 
response plots of filter algorithms created by the PSFSS 
support tool to implement the desired corner frequencies. 
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APPENDIX B 


2920 HARDWARE SUPPORT CIRCUIT 


The nominal specifications of the digital pre- 
sampling filter are: 

1) The filter transfer function is to approximate 
the transfer function of the analog PSF. 


2) 

Balanced 

differential input 

3) 

100k ohms 

minimum input impedance 

4) 

20k ohms 

maximum output impedance 

5) 

Nine bit 

signal resolution 

5) 

Four filters per 2920 chip 


The design specifications for the 2920 PSF filters 
will impact the filter implementation in two basic areas. 
Items 1, 5 and 6 of the specification list are primarily 
determined by the software (within the constraints of the 
2920 chip hardware) the 2920 will execute. The structure, 
the filter algorithm, the accuracy of the multiplications 
and other software techniques determine these parameters. 
The I/O specifications of the filters are determined by 
the hardware. The 2920 will require hardware support to 
meet the input specifications and to perform post 
filtering of the output signal from the DAC’s, Figure 2-1 


C- X 


8 4 



shows a block diagram of the extra support circuits 
required for the PSF application. 

The SCB circuits perform several functions. They 
act as buffers between sensor signals and the 2920 ADC's. 
The circuits convert the differential input signals to 
single line, and act as an anti-aliasing filter and input 
buffer. Finally, the circuit scales the input signal to 
levels acceptable for the 2920. 

The SCB circuits shown in Figure B-l perform these 
functions. The circuit was designed for differential 
signals magnitude limited to ± 5 volts. The corner 
frequency of the circuit is 1/R 2 C. The sample rate of the 
filter is known to be less than 4.166 kHz, therefore the 
corner frequency of the SCB circuit is 2.5 kHz. The gain 
of the SCB's is 0.5 volts/volt. The input impedance of 
the circuit 100k ohms. 

The reconstruction (or post) filters remove the 
high frequency DAC noise prior to resampling by the PCM 
system. The corner frequency of these filters is 6.6k 
ohms . 


The one volt reference voltage should be generated 
locally to the processor to reduce reference noise. 
Figure B-l is the circuit diagram of a circuit which meets 
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the circuit specifications. This circuit was designed 
using available hardware. It is not intended to be a 
space optimal design. The circuit will serve as a model 
for a rough estimate of the space requirements of the 2920 
PSF circuit. Table B-l contains a breakdown of the 
various parts and estimated size. The space requirement 
for the circuit of Figure B-l is estimated to be 10 square 
inches, or 2.58 square inches per filter. 


Table B-l Estimation of the 2920 PSF Circuit Size 
Part Description Size (in^) 


2920-D 

1 ! 
1 Signal Processor 1 

.9 

1 

1 

I 

IC1 

I LM747AH, SCB- op-amp 1 

1 | 

.09 

I 

1 

j 

IC2 

! LM747AH, SCB op-amp I 

.09 

1 

1 

1 

LM110H 

[ Reference Voltage Buffer 1 

.09 

1 

i 

1 

1N746 

1 3.3 volt Zener, Diode 1 

.05 

i 

| 

Crystal 

1 3.2 MHz Crystal I 

1 1 

.57 

1 

1 

1 

1.5k Ad j 

1 1.5 Adjustable potentiometer 1 

1 1 

.12 

1 

\ 

Misc 

I t 

1 24 Misc Resistors @ .12/resistor' I 

2.83 

1 

1 

1 


1 * 
! 17 Misc Capacitors 0 .12/cap ! 

2.04 

1 

i 

1 


Subtotal i 

1 

6.83 

1 

1 

| 


1 

circuit layout factor x .5 ! 

1 

3.42 

1 

1 

1 


Total estimated size (four filters) | 

10.25 

1 


Estimated size per filter 1 

2.58 

1 


APPENDIX C 


2920 PROCESSOR OVERVIEW 

The Intel 2920 digital signal processor is 
specifically designed for signal processing applications. 
The unique architecture of the 2920 gives it the 
capability and the speed necessary to perform real time 
digital signal processing functions. 


2920 Architecture 

As shown in Figure C-l, this unusual stand alone 
microprocessor may be divided into two major sections, 
analog and digital. The architecture of the 2920 allows 
analog and digital functions to be performed 
simultaneously. 

The analog section, which is controlled by the 
digital section, contains all the components necessary for 
analog I/O. This section's circuitry , is composed of four 
analog input channels, an input- channel multiplexer, 
sample and hold circuitry, eight bit A/D and D/A 
converters and eight analog of digital output channels. 
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Data is passed between the digital and analog sections by 
the nine bit DAR (Digital Analog Register). 

The digital section of the processor contains 192 
x 24 bits of EPROM program storage and 40 x 25 bits of RAM 
storage. The 2920 supports a maximum of 192 instructions, 
and up to 40 variables. Arithmetic operations are 
simplified by the capabilities of the shift register, 
which scales the A register operand by 2**N, where: -14 < 
N < 3. If an arithmetic operation causes an ALU overflow, 
the largest valid magnitude of the result and the proper 
sign bit is retained as the resultant. 

The DAR is a nine bit 2's complement binary 
register which passes data between the digital and analog 
sections. After an A/D conversion, the result is 
converted from a l's complement to a 2's complement number 
for the DAR. When the twenty five bit digital section 
reads data from the DAR, the lower sixteen bits of the DAR 
are read as l's. This partially compensates for A/D 
offset errors. 

2920 Instruction Set 

The assembly language used by the 2920 is as 
unusual as the 2920 architecture. Each assembly language 
instruction is comprised of five instruction fields (see 



Table C-l). The contents of these fields specify the ALU 
instruction, the B register, the A register, the shift 
factor and an analog instruction. Table C-l lists the 
2920 instruction set. The execution time of all the 
instructions is exactly equal. 

Table C-l The 2920 Instruction Set 

I (1) I (2) | (3) | (4) I (5) ' | 

I If!! f 

| ALU I B | A | Shift | Analog I 

I instruction | address | address j code 1 instruction ! 


Mne- | 

monies ! Operation 


ALU Instructions 


ADD 

I (A x 2" n) 

+ B > 

B 

SUB 

! 3 - (A X 

2 ~n) > 

B 

LDA 

! (A x 2~n) 

> 

B 

XOR 

I (A x 2~n) 

© B > 

B 

AND 

1 (A x 2~n) 

* B > 

3 

LIM 

| Sign (A) - 

> 

3 


Analog Instructions 


IN ( k) | 
OUT (k) I 
CVTS | 
CVT(k) | 
CNDS | 
CND(k) | 
EOP | 
NOP I 


Sample input channel k 
D/A to output channel k 
Convert sign bit 
Perform A/D on bit k 

Select sign bit for conditional ALU instructions 
Select bit k for conditional ALU instructions 
Reset program counter 
No operation 


Shift Codes 


Rnn | Shift A right nn bits (scale factor = 2~-nn) 
| nn = 0,1, 2, 3, .... 12,13 

i Shift A left nn bits (scale factor = 2~nn) 

| nn = 1 , 2 


Lnn 
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I/O Timing Requirements 

The architecture of the 2920 places several 
restraints on analog I/O instruction sequences. A series 
of IN(k) instruction causes the input multiplexer to 
select channel k for input. The execution time of this 
series of IN(k) instructions must be at least six times 
the RC time constant of the sample and hold circuitry. 
This sampling time allows the hold capacitor to fully 
charge to the value of the input signal. The instruction 
sequence shown in Table C-2 executes an A/D conversion on 
input channel zero. The three IN(0) instructions (3.75E-6 
seconds at 1.25E-6 seconds per instruction) provide the 
settling time required by the sample and hold capacitor. 

The A/D converters require time to settle before 
each CVT(n) instruction. A series of analog NOP's whose 
execution time is no less than 1.2E-6 seconds should 
precede each CVT(n) instruction. The entire A/D 
conversion sequence (Fcry = 3.2 MHz) requires 31.25 E-6 
seconds or 25 instruction cycles (at 1.25E-6 
sec/instruction) . 


Table C-2 Nine Bit A/D Conversion Sequence 


! 


1 

Digital 1 

Analog 

1 

I 

1 

1 

SUB DAR, DAR, R00 , I 

XN0 

1 

1 

2 

! 

* 1 

INO 

1 

1 

3 

1 

★ 1 

NOP 

1 

1 

4 

1 

* 1 

NOP 

1 

1 

5 

1 

* ! 

NOP 

1 

I 

6 

1 

* | 

CVTS 

1 

1 

7 

1 

ADD DAR, KM2 , R00 , | 

CND6 

1 

1 

8 

1 

* 1 

NOP 

! 

1 

9 

1 

★ 1 

CVT7 

1 

1 

10 

1 

★ 1 

NOP 

1 

1 

11 

1 

* 1 

CVT6 

1 

1 

12 

1 

* 1 

NOP 

1 

1 

12 

1 

it | 

CVTS 

1 

I 

13 

1 

it I 

NOP 

1 

1 

14 

1 

* 1 

CVT4 

i 

1 

15 

1 

it 1 

NOP 

i 

i 

16 

1 

* i 

CVT3 

1 

1 

17 

1 

LDA T , T , R00 , 1 

CND4 

1 

1 

13 

1 

★ 1 

CVT2 

1 

1 

19 

1 

LDA T , T , R00 , ! 

CND4 

1 

I 

20 

1 

* 1 

CVT1 

! 

1 

21 

1 

LDA T , T, R00 , | 

CND4 

I 

1 

22 

1 

* i 

CVTO 

1 

1 

23 

1 

LDA INP , DAR , R00 , 1 

* 

i 


— > Indicates available instruction 
— > ^ = 3.2 MHz 

"> C hold - 500 P F 


The sequence of instructions for D/A conversions 
has similar timing restrictions {see Table C-3). Once the 
DAR is loaded with output data, the D/A circuitry requires 
6E-6 seconds settling time before an output instruction 
may be executed. Crosstalk between output channels is 
reduced by executing an OUT(k) instruction, where k is an 
unused channel number. Immediately following this output, 
the signal should be routed by a series of OUT(v} 
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Table C-3 Instruction Sequence for D/A Conversions 


1 


I 


Digital 1 

Analog 

. 1 

1 

1 

! 

LDA 

DAR,OUTP,ROO , i 

NOP 

! 

1 

2 

1 

ADD 

DAR, KP1 , R01 , 1 

NOP 


1 

3 

1 


* | 

NOP 

1 

1 

4 

1 


* 1 

NOP 

1 

1 

5 

1 


* 1 

0UT4 

I 

1 

6 

1 


* 1 

OUTO 

1 

1 

7 

1 


* i 

OUTO 

1 


> Indicates available instruction 

> DC offset correction 

> 

> 


F = 3.2 MHz 

Output channel four not used 



instructions to the appropriate output channel v. The 
execution time of the OUT v instruction series must be no 
less than IE-6 seconds. 

Multiplication in the 2920 

Time delays, additions and multiplications are the 
three basic functions which form the basis for 
implementing digital filters. Time delay and addition 
functions are supported directly by the 2920 assembler. 

Multiplication is expressed as a series of LDA, 
ADD or SUB instructions. In a typical filtering 
application, the algorithms which multiply variables by 
constants comprise 75% of the filter program code. The 
efficiency of these multiplications determines the space 
efficiency of the filter code as a whole. 
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There are several ways to code a multiplication 
between a variable and constant. For example, consider 
the multiplication: 

REGO = REG1 * 1.8727 

where: REGO = variable A in a filtering program, 
REG1 = variable B in a filtering program. 


The first task in the process of coding a 
multiplication algorithm is to express the constant as a 
25 bit binary number. The 2920 has a 25 bit (24 bit plus 
ign) fixed point accumulator; therefore, the binary 
representation of the constant need not be greater than 24 
bits. 


Table C-4 illustrates that there are many ways to 
code a multiplication, some ways more efficient than 
others. Sequence A requires twenty instructions to 
perform the twenty five bit multiplication. Sequence uses 
several coding tricks to shorten the sequence to nine 
instructions. If twenty five bit accuracy is not 
required, the multiplication algorithm can be shortened 
considerably by rounding the constant to an appropriate 
value. See the 2920 Assembly Language Manual for further 
ietails on 2920 programming techniques. 



i 
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Table C-4 Examples of Algorithms Which Multiply REGO 

by 1.8727 


A) Instruction Binary Result 


LDA 

REG1 , 

REGO, 

R01 

1 

1 0.100 

000 

000 

000 

000 

000 

000 

000 

ADD 

REG1 , 

REGO, 

R02 

I 0.110 

000 

000 

000 

000 

000 

000 

000 

ADD 

REG1 , 

REGO, 

R04 

1 0.110 

100 

000 

000 

000 

000 

000 

000 

ADD 

REG1, 

REGO, 

R05 

1 0.110 

110 

000 

000 

000 

000 

000 

000 

ADD 

REG1, 

REGO, 

R06 

1 0.110 

111 

000 

000 

000 

000 

000 

000 

ADD 

REG1 , 

REGO, 

R08 

! o.lio 

111 

010 

000 

000 

000 

000 

000 

ADD 

REG1 , 

REGO, 

R09 

i 0.110 

111 

Oil 

000 

000 

000 

000 

000 

ADD 

REG1, 

REGO, 

Rll 

1 0.110 

111 

Oil 

010 

000 

000 

000 

000 

ADD 

REG1, 

REGO, 

R12 

1 0.110 

111 

Oil 

Oil 

000 

000 

000 

000 

LDA 

REGO, 

REG1 , 

R12 

I 0.000 

000 

000 

000 

110 

111 

Oil 

Oil 

ADD 

REG1 , 

REGO, 

R00 

1 1.000 

000 

000 

000 

110 

111 

Oil 

Oil 

ADD 

REG1 , 

REGO , 

ROl 

1 1.100 

000 

000 

000 

110 

111 

Oil 

Oil 

ADD 

REG1 , 

REGO, 

R02 

1 1.110 

000 

000 

000 

110 

111 

Oil 

Oil 

ADD 

REG1 , 

REGO, 

R04 

1 1.110 

100 

000 

000 

110 

111 

Oil 

Oil 

ADD 

REG1 , 

REGO, 

R05 

1 1.110 

110 

000 

000 

110 

111 

Oil 

Oil 

ADD 

REG1 , 

REGO, 

R06 

1 1.110 

111 

000 

000 

110 

111 

Oil 

Oil 

ADD 

REG1 , 

REGO, 

R07 

1 1.110 

111 

100 

000 

110 

111 

Oil 

Oil 

ADD 

REG1 , 

REGO, 

R08 

1 1.110 

111 

110 

000 

110 

111 

Oil 

Oil 

ADD 

REG1 , 

REGO,. 

RIO 

1 1.110 

111 

110 

100 

110 

111 

Oil 

Oil 

ADD 

REG1 , 

REGO, 

Rll 

! 1.110 
1 

111 

110 

110 

110 

111 

Oil 

Oil 


3> 


LDA 

REG1, 

REGO, 

R09 

1 

1 0.000 

000 

001 

000 

000 

000 

000 

000 

ADD 

REG1 , 

REGO, 

R12 

1 0.000 

000 

001 

001 

000 

000 

000 

000 

SUB 

REG1, 

REGO, 

R02 

1 0.000 

000 

000 

110 

110 

000 

000 

000 

ADD 

REG 1 , 

REGO , 

R05 

1 0.000 

010 

000 

110 

110 

000 

000 

000 

SUB 

REG1, 

REGO, 

R08 

1 0.000 

001 

110 

110 

110 

000 

000 

000 

ADD 

REG1 , 

REGO, 

RIO 

I 0.000 

001 

110 

110 

110 

111 

Oil 

oil* 

ADD 

REG1, 

REGO, 

L01 

110.000 

001 

110 

110 

110 

111 

Oil 

Oil 

SUB 

REG1 , 

REGO, 

R03 

1 1.111 

001 

110 

110 

110 

111 

Oil 

Oil 

SUB 

REG1, 

REGO, 

R05 

1 1.110 
1 

111 

110 

110 

110 

111 

Oil 

Oil 

Note 

: 1.8727 10 = 

1.110 

111 110 

110 

110 

111 

Oil 

Oil 

a. 

2 ■ 

err 


♦ 


May cause accumulator overflow 



Sample Rate Division 


Table C-5 shows a sequence of instructions which 
may be used to modify the sample rate of the filter. The 
variable OSC is a counter which is decremented by the 
constant XPn (n ranges from 0 to 9) each time the code is 
executed. The filter is kept inactive by disabling the 


Table C-5 Sample Rate Division Instruction Sequence 
for Sample Rate Divisor of Three 


Instruction I Comments 


| LDA DAR , OSC , R00 , NO P I DAR = Current Value of OSC 
I * ,CNDS | The time delay and input 

I * , CNDS I instructions execute only 

I * ,CNDS | if DAR (OSC) is negative 

I LDA OSC,KP2,ROS,NOP | Reset OSC 
I SUB OSC, KP1,R05,NOP ! Decrement OSC 


time delay and I/O functions while the counter variable is 
greater than zero. When OSC decrements to a negative 
value then the filter input and output values are updated 
and the time delay function is executed. OSC is reloaded 
with a positive constant KPm and the sequence repeats. 
This sequence effectively divides the original sample rate 
of the filter by the ratio of the positive constant to the 
decremental value. The corner frequency of the filter is 
divided by the same amount. The new corner frequency of a 
filter using the sample rate division may be calculated 


U s ;g the relationship 


c 

F‘ = [ ] 

KPm+1 

INT [ ] 

KPn 


where 


= new sample rate 

F = basic sample rate 
c r 

KPm = positive initial counter value 
KPn = decremental value 



APPENDIX D 


PSFSS 

Pre-Sampling Filter Software Support 

The PSFSS software package reduces the process of 
writing a program which implements four pre-sampling 
filters with the 2920 to a single instruction command. 
The PSFSS command: 

DO FILTER (FI, F2,F3,F4, [$LIST] ) 

links the four filter modules F1,F2,F3,F4 and assembles 
the code. It also produces an optional program listing 
(if specified by the $LIST option) and loads the finished 
quad filter code into the EPROM of the 2920. The chip is 
then ready to be plugged into the hardware for operation. 

The PSFSS package maps the 2920 EPROM program 
memory into four blocks. Each memory block is forty eight 
program steps long. Filter modules are assigned a memory 
block during the linking process. The filter modules are 
written to handle all the necessary filter functions 
within these forty eight program steps, with the exception 
of output. 
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Phase delay due to digital calculations is 
-■nimized by having each filter module output its results 
immediately after the calculations are complete. As 
nodule A finishes its calculations, module B (the next in 
its program space) starts its calculations. While the 
digital section is executing calculations for module B, 
the analog section is executing the output sequence for 
module A. Thus, the output of module FI is handled by F2, 
the output of F2 is done by F3 and so on. 

Digital filtering programs must meet several 
requirements before they can be included in the PSFSS 
library. During the module linking process, PSFSS 
searches each of the modules specified in the DO FILTER 
command for several special characters. These special 
characters aid the linking process by specifying the 
location of various module and position dependent 
instructions. Listed below are these special characters 
and their corresponding PSFSS operation: 

> PSFSS replaces the underline character with 

the block number assigned to the filter 
module used in the code to create unique 
variable names. 

’#' — > The # character is replaced with the block 

number of the previous filter module. Used 
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to identify the variable name of the previous 
modules output. 

‘0UT+' > OUT+ is replaced by OUTk, where k is the 

output channel assigned to the previous 
filter module. 

• IN+ ' > IN + is replaced by IMv , where v is the input 

channel assigned to the current filter 
module. 

* @ * > If the module this appears in is in block 4, 

it is replaced with an SOP instruction to 
indicate end of program, otherwise it is 
removed from the code. 

PSFSS Commands 

The valid PSFSS commands are: 

COPY FILTER. LST TO :nn: 

COPY FILTER. LST copies the assembled list 
file to the line printer (:nn:=:LP:), or 
console (:nn:=:C0:). FILTER. LST contains the 

filter printout and the assembled code 
listing if specified by the.$LIST option. 

COPY HELP.??? TO : n n : 

COPY HELP.??? copies of all the help files to 
the line printer (:nn: = :LP:), or to the 
: = : co : ) . 


console (:nn 
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DO FILTER(fl,f 2,f 3,f 4 [,$LIST] ) 

The DO FILTER command assembles four third 
order Butterworth filters which have corner 
frequencies specified by fl,f2,f3 and f4. If 
the $LIST option is specified, the assembled 
code listing is included in the list file. 
Fl,f2,f3 and f4 are selected from the set of 
filter modules listed by the HELP FILJters 
command . 

DO SIM(FILTER) 

DO SIM(FILTER) sets up the 2920 simulator 
program for simulation of assembled filter 
code. Consult the 2920 Sim ulator User's 
Guide for further information. 


HELP COM 

HELP COMmand lists a short description of the 
common PSFSS commands. 


HELP COP 

HELP COPy lists commands used to copy files 
to the console or line printer. 

HELP EXA 

HELP EXAmple lists a short PSFSS programming 
example . 



103 

help fre 

HELP FREquencies displays a listing of the 
available filter module corner frequencies 
and their respective sampling rates. 

help hel 

HELP HELp lists the available help commands. 

HELP PRO 

HELP PROcedure outlines the proper sequence 
for using the DO FILTER command. 

I 

PSFSS Software 

The PSFSS software is written to operate in an 
Intel ISIS-II operating system environment. The DO FILTER 
command which initiates the creation of a PSF filter 
group, prompts the ISIS-II system program DO (renamed from 
SUBMIT), to process the commands in the command file 
FILTER. CSD (listing D-l). The FILTER. CSD file directs the 
ISIS-II editor to concatenate five files into a temporary 
file for processing. The five files include the 
HEADER.TXT file (Listing D-2) , and any four filter source 
files (Listings D-3 through D-16). The temporary file is 
, edited to link the filter modules, then is assembled by 
the AS2920 assembler, again under the control of the 
| FILTER. CSD command file. When the assembly is complete, 
the program BEEP (Listing D-17) is called to notify the 
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user that assembly is complete. When the user responds, 
the command file directs the loading of the object code 
created by the assembler in the 2920. 

The HELP information is processed by the program 
HELP (listing D-18). HELP dumps the contents of the 
HELP. XXX (where XXX is specified by the HELP XXX command) 
file to the CRT. The contents of the HELP. XXX files are 
shown in Listings D-19 through D25. 
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of POOR 



Listing D-l FILTER. CSD Command File Listing 



FILTER 
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ORlGi?^ E3 
OF POOR 'QU&iUnnf 



Listing D-2 HEADER . TXT File Listing 



ORIGINAL 

OF POOR QUALITY 


*..A ** 4 * 

.'(IF: , 


•00H z 


£*F 


1 101 Oil 


00 110 lie 


uD^ / OUTPf t f\ 0 0 

nLlJ D AnC / ft P 1 / SO 1 

/ 1 u c. r r aCIcNT Cal Cl? la » ion l 

: 31-1 . 47617977c 55 01. Oil 110 Oil 110 Oil Cl 

LDA 31) >11 1 _P 1 * POO 
a j r.i sumi / _p i / n* p i 

51? f> 3UM w _ P i / P.05 
A i> L« Si! ilu.Fl * PO 7 1 0 IJ 7 4 
3U b SUM! /.Pi /Rl 1 sGUT$ 

ci;r» 

'5U t ; DAP » DAP » K*.*0 ? 1*4+ 

; S2-- . o4 1 2S**o65 a ~C*0. 10 1 001 000 010 1*. 

LDfi TEMP ,._T2/P01 / IN* 

ADD 7gHP/_72/R03 
ADD TtMF*/.T2^R06 
ADD 7EMP,_72/Pll 

ADD 7c Mr' * _T2 / r,’ 1 L 

SUB SUM t > TEMP , ROO - £ V 73 

ADD DAp f i\M2 / POO * CND© /ANAcOo 

iAi. AND a2 COEFFICIENTS 

ADD 3iMlf_Zl/R02 
ADD 3UM1> .22 / R03 / CVT7 

.'Cl*. 625540501= 00.101 000 000 Oil Oil 100 010 101 

LI- A :3uf12/-Z?#R01 

AD J 3UM2 / .4.r‘ i FcOo / CO • © 

SUP Sum2>_2P/*i3*CVT5 

NOP 

LOT ^ 

TOP 


. 7 


r;/ 


>57:o-= 00.110 OlO 100 UO 


0 000 00 i 


. c.M ir' ; «Zr* / P'/ i 


• L V 


;.L‘A 

-Dm i j 7 /KOO .• CM 04 
ripD i £Hr'7._c r* .* RO 2 / L : 
v. DA 7*| » KOO / Li\D-4 
ADD i eMP / -IP / P05 / CO 
LDA 7 >i ; r; 0 0 * C M D 4 
ADD TEMP / _ZP # R07 #CV70 

* Pc AD i FOP 1 NPu ♦ 

Li'A INF * DAP / POL 
ADD SUKl/i.NF/RGo 
ADD 3UM2 .• SUtll r R02 
ADD j_cHP / 3oM2 / AO 1 
ADD j_cMP / 3UM2 / ft 02 
ADD • c MF .* SUM 2 * »\ 0 0 
Ai. 5 ! 1 TEMP » 5UJTI t P07 

•TIME DEL A • 

.• *i i" A !'»’ • _ • J i L A y 

l.i * a .DSC / N r 0 / pO "j / C N D 3 
L£*A -i 2 .* -I u SOO / u ?i 1*3 
lDA _Il i INP t POO t CNDS 

L.DA .72 f -P 1 * PO*j £«.N03 

SJP — 03C / *\Pl / R05Q 

hiJS , 

‘.DA uU . P ♦ / * EHP f Lv 1 * CNDb 
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ORIGINAL {*? 
OF POOR QUALITY 


i FILTER _ r c - 15 OH i BP 200 

•_!*« BArC / OUT?# / r!v«) 

AL*0 DAkfKr’i inOl 

# COEFFICIENT CALCULATIONS 

^761/9776= 01. Oil 110 Oil 110 Oil Oil 101 Oil 

LD»i SUM 1 1 J? 1 1 r . 0 0 

»5 d d sum / -.p i * ko i 

o*j D SUM 1 / 1 • POo 

ADD SUM1 » fc r’l » RO? ; OUT h 
SUB Slirtl *„Pt fRlt ,Qd7* 

OUT* 


5 U 5 0 A A / '*4fv # a 00 * I N ♦ 

» B2-~. oh12843o 5= -00. 101 001 OoO 010 IOi LOO 110 110 

.2*4 TEMP * -T2 /POt > IN* 

ADD 7EMP/—72/R03 
ADD TEMP;.T2/R0o 
ADD TEMP r«T2#Rll 
aDD 7 cMP ^ 72 / r: 1 3 
SUB SUM1 /TEMP t ROO » CVTS 

A D D l*m? / /v M2 / t (.-NBo / A iV h L G G 

/A1 pND A 2 CuE* PIC IcMTS 

ADD SUM1/«I1/R02 
Al-D SUM 1 1 .22-/ KOo • CVT7 

• s258*4«)5oi« 00.101 ooo ooo on on 100 010 101 


LDa BUM2/.ZP*R0l 
40 J 3UM2 / .Er'/ r\03 / CVTa 
ADC* SJM2.-.2P/R10 
>'JB SUM'2 / .2P / R 1 3 / C VT3 
NOP 
v T *i 
NC? 


• l / KO 85 . / 9 / 7 25 7 56- 00.110 010 100 110 110 1G0 000 00 j. 

L V a TcliPj; .Zr'i'v'0 1 / CO i 3 
i / T/Rv*/ /CNB^f 
ApD lEMF’ # .Zr » RG2 / COT 2 
wU a 7 .* i •POU/CNOA 
A Dl* TcMF' • .IP * r.O 5 i C V T 1 

•. i * a r » r / poo / cnd a 

D D Vc MP / _ ZP / KO 7 / CVTO 

* : ■> 2 0 r f 0 P. I M PUT 

l.DA I Nr / DAP / *<03 
A D C SUMl/i NP / AO 3 
41*0 SUM'2 / SUttl / R02 
ABB TEMP / SUM 2 / ROi 
A D 0 7 c MP / SUM 2 * K 0 2 
ADO 7 EM P / S UM 2 / 4 05 
ADD TEMP/ 3UM2 / -.0 7 

: “I ME DELAY 


— *A UAi\ * — C S * 'V«v 

i$S 

L Jh --»! / £NP / ROO / CNpS 

L DA . i J / _r 1 » RU*> / CN Bo 

SUB .OSCfKPlfROSG 
LDA JT/ SUM1 / POO .♦ CNOS 
i_DA . Zp 1 3UM2 / ROO .* C ?4 p 5 
-DA *JU • r ' ♦ .• TEMP / LO 1 f uNO S 
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rC=iGOHz 


r «r n * 

c*F 300 


D A pAR t 0!J I r Hr i rvOO 
AL‘0 j AK t Kr* 1 ; K'j 1 

* C£tF r I ClcrJ « CALCULAT IONS 

;B1=1,47*1??776= 01. Oil 110 Oil 110 Oil Oil 101 Oil 

wDA SJrii/.rl/ROO 
ado suni/.Pi /^oi 

SUB 3UM1 / _r 1 / a 05 _ ^ 

AD 0 SUM 1 / _r 1 » R07 / 3 1 J 7 4 
SUB SUM 1 y _r' 1 y ft 1 i y u U 7 0 
OUT* 

SuB D Aft / D Aft / ft 00 / 1 N+ 

\ B 4 — ~ • 6-tl234oo5= “00 • 101 001 000 vlO IvL 100 110 ilv 


v.pA 7£MP / -72 / POi / I N + 

A*JD 7tHr' / fc /KOw 

m 

nDJ i^flr > _T^ t (vl x 

AGO i cMPj _T2 > fti 3 
3U3 SUM 1 / i EMP t ROO i CUTS 

C I 1 GG.\’/ KM2 » ROO »CNDo /mNOLGo 

;Al AMD A2 COEFPrCXENTS 

m or; sum /. 3-i /so 2' 

SUM 1 •/ _22..- R03/CUT7 

;;:i=.ii2534050i= 00. ioi ooo ooo oil Oil too oio 101 

!.DA 3UM2 i _ZP >R01 

a . o ;< sumo _2F /S03 / v r.j 

■AiiU glJM2> _Zr' * RIO 

Su'd SUM2/.*r/Ri3/CVTS 

NOP 

sort 

NOP 


/' ;./.\v=. r‘.-07257So= 00.110 010 100 no no too OOO 001 


:-LiS 


Li/ 1 A 
ADD 


uJ A 

A j u 


T E il P / _ I r‘ t K*j 1 * C U 73 
f y 7 ; ROO y.CN£H 
1EMP#..1P»R02|CV72 
T / r • ROO / CN Cm 
7 clip * «ZP /ft 05 1 CUT 1 
7 / r • POO t CMDa 
* clip f _ZF r R 07 / GOTO 

DY FOR INPUT 


L Da INP ,• B’Aft # R03 
r’t D II oiJMi t INP-/ri03 
ADO SUM2/SUH1/R02 
ADD .cMP/SUM2/K01 
add T£MP-.5U«2^§03 
ADO » EM? / SUM 2 / ROL 
ADD TtflP • SUM 2 $ RO * 

•TIME DELAY 

•_ .j A t i A K > _ 0 ■? * . y K V • ; 


v_ : J A _G‘2>C » f».P2 / kOj / CN U 3 
L D A » Z2 / 1 1 / oO 0 / CN03 
LDA _I1 1 IMP • POO •C?403 
LDA „72 / -Pi / ROO m 7ND3 
SOB .CSC #KPl 9 R05Q ‘ 
uBA »P1 • SUM t yPvO/ONuS 
LDA _ZP^ S‘JM2 / ROO • CMOS 
L if A OU * r TEMP • LO 1 • lN D o 
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ORK5EKAL F-t'Jl 13 
OF POOR QUALITY 


..i 4 i -4 w «. 

H: 


• FILTER - pc«80Hz er 30 

«t ! ‘i i.*n\ / -JjJ I r‘» > KcO 
ADD BAft/KPl *R01 

; COEFFICIENT CALCULATIONS 

• 51* 1.87277499053-#* 01. 110 ill 110 110 ill 000 101 iio 

LOA ?uru /_? 1 /Rlv 
sus oumi .•S um »r83 

A L* 0 3 0 M t / — P 1 # L V t 
SUB SUM 1 1 _r 1 .* PO 3 
3lj S gUM t .» _?1 » n 09 / UU74 
3‘uB bOMl ? _rl /PI 1 / 0UT4 

.! 52 = - . SGo-^iISooOSm 3 -00. lit 000 L01 111 000 010 111 00 i 

LDh TEMP/-T2/K07 »0UT» 

SI/ S OAR t OAK / ROO /IS* 

A 00 TEMP# a T2»R06 r IN* 

ADO TEMP « ?EMP./ SI 0 
SU B T c.MP » _T 2 > Si 2 
ODD TEHP-»_T2/P00 
3JB Tr.MP #™T2 #fi03 
SUB SUMl # TEMP / ftOO / CV TG 

ADO iIIAp. # KM2/KOO # CNt*o » ANALOG 

■A1 AND A2 it Jtr FICIENT3 

ADD SUM1/.Z1/R07 
ADD 3UM1 „• _Z2/SCS#CyT 7 

• Cl = . 38oO?52 < *G?2o= 00.111 000 101 101 Oil 100 100 Oil 

•_SA 3U.M2 > _IP > SCO 
SUB' S'UMU .• _ZP » SOS / C v'To 
ADD SUM? # SUM2 » R 1 3 
ADD SUM? »_Z?> ROB /CUTS 
•SUB SIj’MU t _iP »R08 
SDR SUM2 1 _2P » R 1 1 / CUT •* 

SUB SUMS t _Zr' > R 12 


■ I /KG — . So / 70 oS«l— 00 . I'.'- 

LI* A 7£MP/_ZP#R01 #CV i'3 

^da r » r # sew # cnD“ 

ADD. "cMP.'-ZP/SOS/C'JT? 
uDfl V (T / KOO »CN04 
A;»D lEt*r'. r P ;i n 6 ,cv7i 
wtn _ T l I »RS 5 »CNDh 

i- ; V 

?RS>.liV FOR INPUT 

-l ; m I NP / l*«ft * R02 
hDD 3Utfl/lNP,ft08 
A U I- o Uf12 / S-Jft 1 / ft 0 2 
A Li L» TcttP / 3UN2 / RO 1 
-*jfi T£MP#SiJ.T2*R03 
hOIs i Fi*P>3UN2./R0o 


111 100 000 


oi : in 


* ■ I'cLf* / 

L : ja / .OSC ^ ft 00 

LB A -OSC .• KPO / R05 j CNBS 
L£A _.~2- “1,R00/CN£S 
•»J*i _ l 1 9 1 NP t RvO / LN&b 

k8s :5id;&n8&g M05 

-3A OU TP+ • » cMP f ft 00 / il ND3 
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onmmi page ?s 
OF P0C.1 QUALITY 


*n nn « 

Fc=75H: 


br 300 


L A jAP/OU 7F "9 / f\00 
*-.£ : I* DAR»ii.?i»»\01 

i C JEr r ICitN r CALCSJuAT ions 

: B i =1 . h7o 17?7?6= 01. Oil 110 

uOA 3UM1/-P1/POO 

A l» Li SUM 1 /_r 1 ; f <0 t 


on no on on 101 or; 


oU.Mi / 


.Pi *n05 


ioi ooi ooo oio ioi too no n 


r» z* zi sum » -Pi / so? » out 

I* j ? SUM 1 * _r 1 / ft 1 1 / GUT ? 

OUT* 

V.j t* Z»An* / I* An .• K j 0 / ( N4* 

f 32=- . g **1284365= -00. 

■-If a Tc?5P / -72/?\0i f 1 M+ 

Apt) TEMP / >T 2 / *03 
Vj{V *7£MP*>.LT2/'Ruo 

a 0 3 I Emp / 2 > ftj 3 

•;Uv> oUMl # TcMP ?k*>0 * C'J > S 

•i?D OaR/PMZ /POO mONDo •ANALOG 

* .-i i AN I* A2 C’j£rF I Cl ENTS 

“I * Z ■ 3 UM l / .li.'nOI 
aZ^O SUMi ; ..Z2/P0S /CVT7 

i£ 1 a . 3Z5 c?h'}50.1 " 00 . 10 1 000 000 0 1 1 01 1 

L l l n 3U M2 / -Zr j P.0 1 

?• 0 0 St) M2 t Jlr t & 03 #CVT 6 

Xu'i 3UM2/,ZP/fU0 

“US 3UM2 / _ZP * r 1 3/ CUTS 

NOP 

NOP 


100 010 101 


». / s 

790725756* 

00 .110 

r. - 

TEMP t _ZP / PO 1 / 

Wi: 

— A 

f / 7 » KOO/CND** 


nM 

TEMP ; — Zr* / 7<02 > 

CVT 2 

lI*A 

• / r y kco / cn c h 


A lip 

TEMP # _.ZF’ r P05 .* 

COTi 

Li* A 

It T y r! 0 0 / C N D 4 


mOu 

TEMP / -Zr t POP / C T o 

r. X A 

nr fop input 


Li i»-. 

inp t Oar / r<os 


r-I'j 

SUMI .» INP/POS 



c*i JM'2 * SUM 1 / k 0 2 


.< i.l |.l 

TEMP t SO M2 / RO 1 


A Li L' 

TEMP 7 SUM2 .» P02 


A J ‘J 

:lMP7SUM2/R05 


“li. 1 L* 

TEMP t SUJ12 7 P 07 


■' .n; 

1 -.'tiLA \ 


L . • s 

Jnh! • -OoL *nOO 



too 000 0 01 


cJh _0b C / KP3 • P*;li • *_ N L>S 
LI- A -22 / -Z 1 /ftOO/CNDS 
—21 f lNP/nOO#CNDS 
L 1*4 _ f 2 / _r 1 / POO / C Nil r* 
oJP _ J3C / KP i / ».'5y 
L I* A - Pi -• SUM 1 .• n 0 0 • C N Z‘3 
LZ»A — Zr » SUM 2 # POO » CiiPS 

A 3U ; ?+ » TEMP , LO i .• CN03 
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I 44 ; 4 444* Y *.*»f;*.*N*.4444*4*r4.4 4*44 4 444444*444w*4^-*44444 4*4-**. 

» •“ iL i c.R _ *C s 70Hl BF 7*) 

LDA BAR /0U7P9 * &00 
ADO JAR * KPi » #40 i 

.• COCr r IC ItHT calculations 

:dl«l.i3l>?356203» 01. Ill 000 111 010 101 110 001 OiO 

wDri iUMl /^.Pl/ftOo 
31: 5 SUM 1 / _.r’ 1 / ROy 

D D SUM 1 # _r’ 1 / R 1 2 
A j D '3 U M 1 * SUM 1 1 R05 
ABB SUM 1 / —PI pi-0 i 9 GUT 4 
S'./t* SUM 1 / ..r’t /.R03 f *JU * ^ 

;B2*-.89V«Ui52« -00. Ill 001 100 110 000 010 010 100 

Ll'A r£KP,_T2jR05/0UTJ> 

SUB BAR* / BAR r ROQ / IN+ 

'ji)B 7E*1P * ~T2fR07 r IN* 

ADD TEMP > TEMP »RG4 
AuD 7EMP»-72/R0v 
SUB TEMP i^T2/P03 

MOP 

SUB SUMl fTEKPrROO# CUTS 

AdB l!AK’ / N. i2/^00 .* CNu-o /ANAuOlS 

•A1 AND A 2 COEFFICIENTS 

ADD .SUM1 » -21 /R07 
ADD SUMl > _Z2 / R08 r CVT7 

;i;i = .svyoA5Bo7= 00.111 001 100 100 in 100 no 001 

LDA 3UM2 /_ZP t ROO 
SUs bUfl2 » — i.F‘ / ROA / L 
ADO 3Utl2 / 3UM2 / Rl 2 
SUB SUM2/-2F /R04/CVT5 
hDl» SUM 2.- -IB *RC'^ 
a 0 D S JM2 / _ £ P / RO 7 / C V* 1 •* 

ADD 3UH2 1 -IP /RIO 

n/K0=.36?3o?o-f9 = 00.110 ill 101 000 111 100 000 010 

LDh i EMP / _Ir* / ROO t CUTS 
..DA T » 7 -ROO/CND* 

SUB TEMP / -.ZP / RG3 / CUT 2 
a Da • * T * ROO / CND A 
SUB TEMP / _ZP / ROS / C V ^ 1 
L DA i • T f ROO / C N 0 A 
CU 70 

•READY FOR INPUT 

EDA INP/DAR/R03 
ADD SUMl / INF* *R0:$ 
a i* i.i 3UM’2 .• SUMl t R 0 2 
r D D *_ E M k / 3 U M 2 / R 0 0 
- i i ~ MP / SUM2 / R 03 

•: J : c. MP t SUM 2 / r. vd 


_ Da JAR .• _ 23C t KvO 

L J A _03C t Kr 0 1 ROS / UNDS 
fc-DA « i.2 .* ~a 1 / ROO / u NuS 
LDA .Zl#lNP#R00tCNB5 
IDA _72/-Pl*R00 -CNOS 
SUB -OSC f K?1 .* fcOSS 
L DA _£* / JL'MJL /. R v 0 1 0 N 0 3 
i-JA .ir’ / bb?1.. / / lNl ; 5 

LOA OUTP + • TEMP * ROO * CNiiS 
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ORfGcM. PAGE' TS 

OF POOR QUALfTY 


tfr 


Li'n Jr>n *£MJ7r£ / ROO 
*00 

* COEr I TIONo 

\ B L = 1 . -iTol/^'To 52 01. Oil 110 Oil 110 01 i Oil 101 Oil 

-0* Surtl/.Pl^gOO 
>iub aUm /^rt ?rCOl 
j U3 SU?1 1 .• Jr i / 03 
*1*0 SUM! /../I ja07 /OUT -♦ 

SUB aUMlf .Pi/^ii/OUT# 

OUT 4= 

SUB OAR ? G*r£ rRGO / IN*’ 

;BH=-. o*lSil*365e -00. 101 001 000 01 u 101 ICO :i-> lie 

.Ll A T£MF*?.,T2# R01 / IN+ 

ADO TZMP>..7£>f;03 
ADD- T cMP *..7 2 / ROo 
ADD TEMPf-J^Rll 
ADD TE«P/_T2/R13 
sub sum i ? temp » roo . cvrs 

ADD OAR /^MZ/ROO / CNDo fiiNnuOu 

/ 6i AND 62 COEFFICIENTS 

A 0 D SUMl/ ...3 1 1 R 0 2 
ADD - SUM 1 »U?2 / R03 / CO T7 

/i:: = .o253AU501= 00.1)1 000 000 >11 OH 100 010 101 

£ 

LDA SUH2/.ZP/R0 1 

ADD' SUM'2 /_£?•» R03 / C'J T o 

ADD SUM2? _.ZP /RIO 

SUP 3UM2#„ZF‘/Rl3/CVTf5 

MOP 

i:u in 

MOP 

••1/Kya. 790725753= 00.110 010 < 00 110 i. 

I..DA TEMP/ _ZP/ KOI ;CVT3 
u.<A 1 » i /ROO/CND*/ 

ADO TcM?._ZP.F;02/CVT2 
i_DA T« if ROO/CNDa 
AD 0< TEtSPrf ,lZP /R05/CVTI 
LCA r,r#R50/GN04 
ADD TEMP»_ZP»R07.CVT0 

; RE A or FQn INPUT 

LDA IMP »DAR t RG3 
ADD SUMl / IWP/R03 
ADD SUM2. • SUMl /KOI 
ADD TEMP /SUM2/ KOI 
A OD TEMP / SUMS / ft 02 
ADO TEMP / SUM:? /R‘)5 
ADO TEMP/3UM2/R07 

»T?ME iiEuA* 

.. D A i.i A ft ■ _.0SC # F. 0 0 


L 0 j 0 0 0 0 0 !. 


lL'H 
u DA 


.0SC/KP4 / R05 / CiUi'o 
‘?»_2l /ROO /ENDS 


LDA _Z1 /INP/ROO »CND3 
c DA _ ; 2 / _?1 /ROO >uMD3 
SUB _0SC/KP1 /R053 
EDA .PI/SUM1/R00/CN0S 
LDA _ZP /SUM2 • ROO /CNPS 
EDA OU TP* / TEMP -LOl/ CMOS 
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• FILTER _ 

ul'A *. | r i K# OUTPt *r<0') 
iiD & P Arv 9 i\P 1 * KO i 

; f.OZF - l C 1 1 N T CALCULATE Of IS 

/ 51 =1.4/61 7? 776= 01. Oil 110 Oil 110 Oil Oil 101 

t. 0 A SUM I. / -? 1 * A 00 
ADD cUMl / «P 1 /R01 
SlJ 5 SUM 1 / _•*’ 1 / i\'J5 
Ap 1* Suftl .• .P 1 / 50.7 / uu • - 
SUB 5iJMt/_P!./f<ll;uUT^ 

GUT* 

•2 U t* 0 A K » i* A ft 9 R 0 0 / 1 N ♦ 

• t‘2 - ■ . 6 h 1284365 s -00 . 101 001 OCO 010 101 100 110 

ui ; A rcMP/.T.2^0l’/IN+ 

.-i DO .TEMP* _?2<gRG , w 
ADO TEMP /-T* >*06 
ADD i EMP/ .TtiPl 1 
ADD TEMP/O.T!!. Pr<1.3 
SUB 3U«1 YTEMP / ROO • CVTS 

A OD J A p,'K^27 ROO * C N Go /ANALOG 

f A l AND." A2 tCdEFFIG IENTS 


Oil 


i !•: 


ADD.: SUttiy.LZl *R02, 

ADD sum /^Z2> R03 / CU t 7 


;G: = U25840501= 00.101 0*00 000 Oil Oil 100 010 101 

L DA 3UM2 f -ZP t RO 1 

a D P SU N2 1 _ 2 P / ft 03 / C V To 

ADD 3UM2/-IP/R10 

SUB SUM2 / -ZP / ft 1 3 / CUTS 

NCP 

CVT-# 

MOP 


H/.<0=. 790722756= 00. HO 010 100 110 

LDA TEMP / — Zf* * PC 1 .* CVT3 
... da . .* r .■ K >0 t CNlH 
ADD TcMP / —ZP t R02 / CL* 72 
L D A » / T .» K 0 0 / C N 0 4 
ADD T EMP ? -ZP t RO 5 > CUT 1 
L i* A i / » »R0 0 / L N D *♦ 

ADD TEMP t -ZP / R07 / CUTO 

• READY FOR INPUT 

L Da IMP r PAR • R03 
ADD SUM1 .» INP *AG3 

M fgSp;lt3a-l;|Sf 

ADD f^MP / 3UM2/ RC’2 
ADD TuMP/3uM 2/R03 
-DP TEMP / 3L'tt2 / RO 7 

• TIME DELAt* 

i. A * ♦ _ 0 S C .• R 0 

L D A _ 03 3 / * P5 / RO 5 / v MI* *? 
uUA — — 2 .* 1 /ROO/ uNl'-a 
LL 1 A — Z 1 / IMP / ROO/CNCS 
LDA -T2/-P1 /R00/CND3 
3*J £> — OSC t KP 1 t R05& 

LDA -Pl/SUMl/ROO/CNDS 
LDA — ZP / GUM2 9 ROO • CND3 _ 

LD A OU • P + / f EflP / *.0 1 / Crl i- 1 3 


100 000 001 
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PAGINAL i> tfle . . 

* Po °» 


: FX1.7ER -. re -'-40 Hi Br o’* 

EDA i*AP / QUTf'S a r!00 
A D D BAR t l\f’X i «v0 1 

,'LJirr «C^Crii LiiLCJLm • -i UNo 

;£l = i.87277-*990534 = 01. 110 111 110 110 ill 000 LOl 

-BA 5 L'Hl / _? 1 / R 1 2 
SU ^ SUM 1 SUM! .► k 03 
ADD SUM 1 / _Pl .»c.0 1 
BUB £UMl.-«Pl*r<08 
SUB 3UM:/^Pi/^:09/0U74 
d'u o SUrti / — P 1 ; RI 1 / CU » ♦ 


• y» = “« 3Sb4359oS5*> s — 00 • 
EDA TEMP ? _„T2 -• E07 / OUT $ 


111 >00 1 ■> 1 ill 000 0 1 ■-* 111 


■ju 3 GAP » UA& » HOO 1 T N ♦* 

ADD TEMP>.-T2,R0S/IN+ 

ADD TEMP > TEMP >Pi0 
SUB TEMP *_T2iR12 
ADD TEMP y~T2.» POO 
SUB TEMr\-„T2 * R03 
SUB SUM!# TcKP.-t\00/CUT3 

ADD BAR t KM2?R0*‘) » CNDo t ANALOG 

/At AN D «c CQEFr i Clears 

rpp SUM! .•-.71 #R07 
ADD SUM 1 / / f\08 / l V f7 

;C: = .83o095240726= 00. ill 000 101 101 Oil 100 ICO OS 

CD A :3L»M'2/„ZP/KCO 
SUB SUM 2 ✓ -IF" r GO 3 / CUT o 
ADD SUM'Z / SUH2 » R 1 3 
ADD SUM2^ZP*R0o/ CUTS 
SJB 3UM2 / -ZP r ROB 
BUB SUK2 f _ZP > R 1 1 r COT B 
BUR bUM2 / ..Zr‘ .• n 1 2 

* l/VvU— . 63 S’ i © 77G352- 00.101 000 ill 100 000 ill Oil ii 

CD A TEj^P * _IP / RO 1 / CUTS 
CD A T • T / RGQ / CND4 
AjD * EMPj— zP .• R03 / CU 72 
!_ y A T * T • k 0*> * uNiA __ 

ADD TEMF’ * —ZP • ROd / CUT 1 
CD A i t T / ROO / CM J-4 

;-V? ro 

* READY FOR INPUT 

..DA INP .• DAR/R02 . 

ADD BUM l / 1 Nr r ti-J cj 
A 0 If SUM2 / SUM 1 / RO 2 
ADD TEMP / 3UN2 .• ROi 
ADD 7cMP/SUf12/R03 
ADD 7EMF* / 3UM2 > ROc 


• • • > >z. fuun \ 

— L A tints t _0oL .* A'J*.' 

COa —QSCjKPI /R 0t> »CNpS 
EDA / -»Z t / R00 .* CNuS 
CD A _Z< tZNPfgOOfCNQS 
C DA «. i 3/ „P1 /cC90/CNJ)3 
SUB -03C.*KP1 #R*>5»> 

*^DA -PI / 3UM1 /ROO/CNDS 
CD A —Zr t SUM2 f POO t Cm Do 
DA OUT r •*■ t i cMF / KOO/CND: 


Vifi 

t~.Z 


ii-;- 


091 
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; FILTER _ r c =35Hz of* ?•;• 

LDA DAni / OUTPfc rROO 
hOU IlAfv^KPl/tvOi 

•COEFFICIENT CnuCUL*riCNS 

#Bl*i. 889336203- 01. Hi OCO 111 01) 101 110 OCi 010 

LDA SUM 1 r J? 1 / rfOo 
SUB 3UM1 / _P1 /r’yi 

* B i* SUM 1 / _P 1 ;ri!i2 

ADD SUM 1 t 3UN1 / R05 
ajCi SUM l t^rX / L.0 1 .• 0 U 7 -* 

SUB SUM! #_P.l/ft03/0UT$ 

;dO*-. 8?99 U152* -00. ill 001 100 110 000 010 010 100 


LL>4 TEMPi-72/R05*0UT* 
oi;S BAR / OAR r AOO / IN* 


2 * RQ7 / IN-*- 
cn.c-.2To4 


sve : IEE3S-' 

njj . i Cl ir / I wl 

ADD TSMPi-Ig/RQQ 
SUB- r£MP,_t*»ROj 
NOP 

StjB SUfll > TEMP/POC * il‘- | :; 


ADD OAR / KN2 .*RCO / CND<6 {ANALOG 

.• A 1 AND A2 COEFFICIENTS 

ADD SUM1 R07 
A'.‘D SUMl / _i.2/P08 > Q'J ~7 


: C 1 - . 89? 645867= 00.111 00 1 100 100 111 100 110 001 


SLID 3UM2 # _ZP > Ro2 » CVT5 
ADD 3UK2.3UM2/P12 
SUB SUM2 j_ZP»R0a .C0T5 
4 u [* $UM2 / -ZF/ROo 
ADD Su«2»_iP/R0/ .•COT-. 
rnJv iiUMZ • K 1 0 

f 1/KO— . SoSSo'r g 4S‘ = 00. i i 

uDA TE*1P/_2P.»P00 mCVT3 
‘ 7? T/R00/CND4 

TE«P/_ZF’,a03/CVT2 
T^T » RGO / CND4 _ _ 
TEMP / _Zr‘ f ROo / CvT 1 
7/7 /R*OG t CNB4 


LDA 
SUB 
I.Dn 
SUB 
LD"t 
£V 7 " 


ill 10 1 00 •) lit 


•RErtt'Y FOR* INPUT 

•_ L‘ 4 .i Nr .• D r»R / RC »} _ 

ADD SUM ! 1 INr / R06 
A l» D SUMS / SUM 1 / R02 
ADD TEMP • SUM 2 1 ROO 

lIs?'3MCS^y3 

JUD Icnr / 

y'l'f-XL DELA 

L i* A ’-A." * _ 0 SC KO0 

l-BA .OSC/KPl /R05/CNBS 
LI ; h JZ2t _Z1 • ROO /CMOS 
LD« _Z1 .• INFsROO/ImJvS 
LDA .72/ -Pi ?ROO /CMD3 
SUB —CSC # KP 1 .• R05y 
L04 _? 1 / SUM! .• ROO / CNOS 
LDA JZf-’ > SUM2 .• ROO / CNDS 
LDA OJTr+j » cil* /RvO/CNBS 


ICO ‘'»00 010 
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1 


FILTER _ 


F c=20Hi 


:iAhf.yj7Pt‘ROO 
ADD DAR>KPl/r\Ol 


•COEFFICIENT CALCULATIONS 

SB! =1.372774990534= 01. 110 111 Ho 110 1 

LDA SUfU/_Pl,cl2 
SUB SUM1/SUH1/R03 
ADD SUM! 

HUB 3UM1 / _P1* R03 
sya .SUN! t _P1 * R09 .• OUT -t 
SUB'SJMl ._P1»R11,0UT* 


38343563854= -00. Ill 
LDA TEMP * _T2 > R07 » OUTS 


000 10 l til 


3U0 ' OAR jDAR »R00/ IN+ 

ADD TEMPO -T2 i R08 I N ♦ 

ADD TEMP VJEMP* RIO 
SUB TEMP'/LT2j R1 2 
ADD .TgflPO._T2>POO 
SUB •.\T*Mp>l2,iS^R0.3 
SUB SUM! 5. JcMRV aoovev TS 

ADO* DAR#KM2VR00 VCMOo i ANALOG 

• A1 . AMO :«2.#okfFICIENTS 

ad li IumIvlM J IS e / c v r 7 

SCI*. 36609524072 <>= 00.111 000 101 101 Oil 

LDA 3UM2> _ZP»ROO 
SUB SUM2,_ZP,R03iCYT6 
-■•DO. SUM 2 > SUM2 / ft 1 3 
ADD SUM2# Jl?> ROo , C0T5 
3U3' 3UM2/LZP * ROij 
SUB eSUH2#*_ZP»Rll /CUT4 
SUB 3UM2 »_ZP-» R 1 2 

J l/KO= . o3?o770352= 00.101 000 ill 100 000 

LDA . TEMP > _ZP / RO 1 / CVT3 
LDA i > • •ROO»CNDa 
ADO TEMP f-ZP/ROS/CV T2 
LDA T >'T > R00 »CNB4 
ADD 7EMP/_ZP/R0o»C , .'Tl 
LDA 7 > T »R00 f CND4 
CVTO 

STEADY FOR INPUT 

LDA INP,DAR/R02 
ADD SUM1 t INP i R08 
ADD 3UM2/SUM1/R02 
ADD TEMP ;SUM2fR01 
ADD TEMP /SUM2/R03 
■•DU TEMP .• SUM2 1 ROo 


•IME DELAY 
-DA DAR .< _C3C j 


too 


LDA _CSC »KP3»R05/uNDi. 
LDA _Z2/'_Z1/R00/CNDS 
LDA _Zl>INP»ROOfCNDS 

bee 

h& 

l.P A QIS it**,, uMr , ROO # 


r»~ o 'j r.\. 


i: oo o id no 


00 010 ill 00 L 


100 100 Oil 
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ORiGS&Al PAGE IS 
OF POOR QUALITY 


rC = ioHz b r 


* ~ * * * 

Mr 


• 4 4 .♦ > ** "* ■» .4 w 4 si '<• 

* “li-TcR 

i^DA DA^/OuTr'I/ROO 
ADD DAR / KF' 1 1 RO 1 

/COEFFICIENT CALCULATIONS 

/ 31=1.37277*990534= 01. 110 111 110 110 111 000 L01 110 

LDA 3UM1 #«f'l /R12 
aii? SUffl /SUN.l / a*</3 
A J J aljfl j .r _r‘ 1 f uOi 

5 J3 Slifll / -Pt /R03 
3U3 -S?JHl#_Pl#ft09f0UT~ 

r)U 3 . 3y«l / «P1 7 R 1 1 / OU T if 

/ S2 = * . 83'6‘h856o 85^* _ 00. Ill 000 

v.DA 7S«P#-T2#‘R07iOUT# 


ill , A'V J- x - v 


SUB DAR t BAR # ROO / IN* 

ADD 7D1P,_T2/R08/IN+ 

ADD 7EHP/TEHF >R10 
SU3. TErtPV- T-2SR12 

‘SUB; sum /temp > roo , cots 

A jC yAt< > «rt2/R00 »CNOo ! AfMlJS 

Jftl'MJD A2 COEFFICIENTS 

. ADD' ,<_Z l/RO? 

ADD 'aU*il / R08 >i»V i 7 

;Cl = .S3o0?52-«072 6 * 00. ill 000 101 101 Oil 100 100 )U 

LDA SCM2 * _ZF » ROO 
cUS SUK2 / _2P »R03 #CVTo 
ADD 3UH2/SUM2/F:13 
ADD SuN2 / _2F / <"<0o *C ! »' f .0 
SUD S<JM2 1 _xP»R08 
fUP SUM 2 / _2P / R1 1 /C’OT-* 

S'dD CijM2 > _*.P .• n'12 

l/;<0=.ti3?c770352= 00.101 000 111. 100 000 ill Oil ill 
Lia TEMP »_2F/R01 #CV73 


'-f? 
-i uEi 


UT,ROQ-CUai 
TcKP f _ZP/R03>CyT2 
-DA T* T trvOOfCNI'4 
nilti T w NP t _IF' t KObtlOll 
Li- A 7/ T .•ROO/CND-** 

CVTO 

fftZiiDY FOR aNFUT 

LUa IN? / DflR / rlC'2 
ADD SUfll/IVP/ROS 
ADD 5'JM2 i SUiil i R02 
Apt! 7£^p r SUM'S i SO i 
AuD TEMP / 3»JM2 1 ROo 
ADD . EM?' t 3UM2 • ROo 

• TIME Dcaa"; 

..DA CAR » _GSC / ROO 

EDA _QSC»KP4.k 05 -CfUj'S 
_DA _Z;;,_il/ftod/CNaS 
LD A «ii iINP/ROOiCMIS 
LDA .72/J>l#RCv.'CNDS 
303 — OSC /Kr'l / R050 
ll'A .Pl/SUrtl/ROOfCNDS 
l. 3A —Zf .• SUrlC .* ROO # CtJDS 

Li* A QU7P+ r TEMP *R00 .» CNDS 
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1 IS 


j 

i 


* rXi-Tcu _ Fc=14Hl or /*■ 

'_ pH phft / .qu Tr <> t ft oo 

Ajj Drtft / Kr'l / r<0 1 
! COEFFICIENT CALCULATIONS 

;Bi=i .G3?33c202= 01. Ill 000 111 010 101 110 001 010 

Li* A SUM 1 f ~P 1 / RO<S 
SUB SUM1 >_P1,R0? 

ABu SUM l / -Pi / R12 
ADD SUM1 >3UM1 ,R05 
ADD 3Urtl/-Pl/L01/0UT4 
SUB SUM! > ..PI / R03/ OUTft 

» 32=- . 3*?9'/t 1 152= -00. Ill 001 100 110 000 010 010 LOO 
LDA TcMP»_T2/P05/0UT* 
oud DAR*DAR#ROO/ IN* 


nSf iSMg j _I 2 > RQQ 
blit! I £:'*P / _ i *. / r<0 J 


nl'k' PAR / ftM2 .• ROO > CNDo i ANALub 

> A 1 AND A2 COEFFICIENTS 

ADD 3UI11 »_Z1 »R07 
ROD 3iJrtl._22;R0a»CVT7 

;<;i=.399645So 7= 00.111 001 1 oo 100 ill IOC 110 001 

LDA i/ul!12» _ZP /ftOO 
SUB SUM2 ; _ZF <R0** >C0T6 
rl'Jj SUM2/.3UM2/R12 
SUP 5UM2 r _ZP »R04 > C0T5 
syp SUM2 / /ROg , 

AiJU SUM2 t _F‘ i RQ 7 t cv r * 

SOD SIJM2 / -ZP » RIO 

• i/l'0= . ■'io?269os3= 00 .110 111 101 000 111 100 000 OlO 

'_PA TEMF' / _ZP /ROO / CO ”3 
'_£lA T/T/P:Q0/CNI)4 
SU6 TEMP .< r' / R03 / C V 72 
t-DA T > T /ROO /CND4 
dUB i'EMP / _Zf‘ / R03 / CVT 1 
lIs 7 / T / ROO » CND4 
CV TO 


/READY FOR INPUT 


LD<S 

hDD 

A D »J 


INP,DAR/R03 
3UM1 / IMP f R03 
SyrtG/SUMJL »k-02 
i clip r SUM2 1 ROO 


3 <.'J F mP f q LJM 2 / K 0 J 
j • J c- 7 c Mf* .* SUML t K* : j o 

• ! 1 Me De-A i 

wi-A J A A f _OSC / R 00 

LL‘A -GSC / KP4 / R05 / CND j 
uDh _I2> _ZL /KGO/CnDS 
LDA _Zl / INP/ROO/CNOS 
- T2 t _?1 9 ROO/ CNDS 
SUB -0SC#KP1 #R053 
u3A ,Pi/SUMl/FOO/CNpS 

LB A — tr • SUM'C f RoO i LNjjS 

L DA 0U7P*/ 7EMP /ROO /CNDS 
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ORIGINAL PA£E !3 
OF POOR QUALITY 


• rxLTEF: _ Fc=10riz cF 60 Hi 

LPA ?iin , *'jjjTrft 1 '^0') 

ABB jnRlKrt/^01 

; coefficient calculations 

» E>i-i . 37277-*'? < ? I )53 i = Ox. 110 111 110 110 ill 000 101 110 

LCA 3UMI/-P1/S12 
SUB SUM 1 / SUM 1 / R03 
ADD Sijflt i _P1 itOl 
•:iLjl» SUM! / _F’l / R03 
iJS SUM 1 / _rl t POP 1 0UT4 
SUB SUM 1 / _Fl » SI 1 » OUT* 

;SS=-.36o435o<i85-»* -00. ill 000 101 111 000 010 111 001 
LBA TFM”‘-_T2.-r!07 »OUT» 


SUB BAR.iDAR>A'00/XN* 

ADD IEMP<_T2»R08, IN* 

ADD i c“P < TEMP. RIO 

SUB 'EM-' • _T2»Rlx 

ADO 'EMP / .72 1 ACO 

3'jS TEMP * _T2 j P03 

SUB Si-Ml t •’Ct;-'/ RO0/CV73 

ADD Da<* .• fi.M2 / ROO mINDo /ANALOG 

;ai and a2 coefficients 

Avi' 3UM1,_S1/P0J 
ADS 5iJMl/_Z2»S06/CU77 

•Cl*. 386093240 '26= 00.111 000 101 101 Oil 100 K 

5«JM2 / -ZF * ROO 
SUo SUM2>-ZR.,R03/CUT6 
A D D S • j M2 1 3 Lfl 2 1 R 1 3 
ADD 3Utt2p-ZP/g0p#CVT5 
SUB 3UM2/_^P/R C*d 
3 US 3UM2 > «Zf' / R1 1 / CVTs 
•SUB SlJrt2# _2P>R 12 


J I/K0*.o39e?70352* 00.1 

LI- A 7EMF* .» _Zr / RO 1 .» CVT3 
LI* A 7 r T /RGO/CNBh 
ADD 7 EMf* _ ZF* / RO 3 / C V T 2 
LDh 7/7 /R00/CND4 
ADD TEMP / -IP / RC*6 # CVT 1 
;j{^ 7 / T .• ROO f CND4 

JREAD'f FOR INPUT 


ooo n: too ooo ;ii on m 


: r gA 
A Is D 
ADD 
Ai« D 
ADD 
ADD 


I Nr • DA R / R02 
SUM 1 / I Nr / R03 
SUM 2 t SUM l / RC2 
*7 cMP t SUM2 / R‘J* 1 
TEMP • SUM 2 / R’03 
7EMP » SUM 2 1 RC-o 


/TIME Dc.i_.Ai* 

•_ L* A j«R » _Ob'v # K'y \j 

LDs _QSC * |>F' 7 » R05 / C NT'S 
\ _i*A — J * — Z1 / ROO / CMOS 
LDA _I1 / IMF *R00 1 C«4D5 
LDA _ 7 2 / _P 1 » R 0 0 * CN 0 3 
SUB — OSC .» KP1 .* A’OjI* 

Ll*A -F 1 * SUM 1 * ROO# CNi)5 
LDA ZP; / SUM 2 • ROO , QN^S 
i .DA OU P-r# ifiM^/KOv/CMOS 
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LINE SOURCE STATEMENT 


1 s 

TITLErfiLARrt PROGRfiN 

2 

ktQ^C 

ftnm 

BEEP 

3 

4 ; 

EXTRH 

CO.CSTS. ISIS 

5 

CSEG 


6 BELL 

EQU 

007H 

7 EXIT 

8 ; 

9 START: 

EQU 

9 

18 

CALL 

CSTS 

ii 

RRC 


12 

JC 

DOW 

13 

HVI 

a BELL 

14 
13 ; 

16 BEEP: 

LXI 

H0C980H 

17 

IW 

L 

18 

JNZ 

BEE? 

13 

cai 

CO 

28 

iw 

H 

21 
22 ; 

23 NBEEP : 

JNZ 

BEEP 

24 

25 

26 LOOP: 

LXI 

HL60H 

27 

IW 

L 

28 

JNZ 

LOOP 

29 

IW 

H 

33 

JNZ 

LOOP 

31 

32 ; 

33 DONE: 

JW 

START 

34 

HVI 

C/EXIT 

35 

LXI 

D.EBLK 

36 

CALL 

ISIS 

37 
23 ; 
39 ; 

JW 

START 

40 

41 ; 

42 ELK: 

DSEG 


43 

44 ESTflT: 

DN 

ESTflT 

45 

46 ; 

OS 

2 

47 

EM) 

START 


Listing D-17 BEEP Program Source Code 


ORIGINAL 

OF POOR QUALITY 


UK 


sower 

sTATOfEWT 

UK 


SOURCE 

STATDtOiT 






HAAS 1 






75 QSJt: 



l 


<ELP 

55 


3X2 

EXCEPT 

74 

08 

QAET 

2 


EXTE) 

CO/ ISIS 

54 

WITF: 



77 

88 

of Hi 

3 




57 


UI 

H/IIUE 

78 ACCESS: 

88 

i 

4 


CSES 


58 

loop; 



77 ECHO: 

88 

9 

5 




5? 


80V 

C/8 

109 

88 

OSTAT 

4 

EXIT 

m 

7 

40 


CAU 

CO ' 

191 0AET: 

SS 

2 

7 

KM 

Effl) 

3 

41 


KV2 

A/0EEX 

102 OSTAT: 

es 

2 

8 

OfO( 

EQU 

9 

42 


CXP 

8 

193 CEILE: 



7 

CLOSE 

EfflJ 

1 

43 


22 

CLOSE 

104 

89 

’:F9:HEL?. ' 

19 




44 


Itt 

8 

105 EXT: 

38 

'KL' 

U 

REABC: 



45 


3P 

LOOP 

194 C8SJT: 

es 

14 

12 


AVI 

C/REAS 

44 

CLOSE: 



197 

88 

OfEH 

13 


UI 

3/CSU 

47 


IWI 

C/CLOSE 

198 



14 


CAL 

ISIS 

48 


UI 

8/CL2LX 

197 RBU: 



IS 


lba 

8STAI 

47 


CAU 

ISIS 

119 RAFT: 

BS 

A 

c 

14 


m 

A 

70 


Utt 

CSTAT 

111 

88 

I8UF 

17 


M 

EXCEPT 

7L 


m 

A 

112 ROT: 

88 

2041 

18 


UI 

H/CSWF 

72 


M 

EXCEPT 

113 

88 

actual 

IV 

SPU: 



73 

mx£: 



114 . 

88 

RSTAT 

29 


WN 

A<8 

74 


891 

c.exit 

115 ACTUAL: 

3S 

2 

21 


CPI 

088 

75 


UI 

I/E8U 

U4RSTAT: 

JS 

2 



22 

(POF 

74 


CAU 

ISIS 

117 IJUE: 

SS 

2048 

5 


CPI 

20H 

77 


JXP 

OPOIE 

118 CJU: 



24 


282 

WE 

78 




117 

88 

1 

25 


INI 

8 

79 

EXCEPT: 



129 

88 

CSUE 

24 


JK 

s?a: 

89 


U! 

H/8UE 

121 

8V 

14 

27 

WE: 



81 

UJP: 



122 

88 

ACTUAL 

28 


Itt 

8 

32 


80V 

C/8 

123 

88 

RSTAT 

27 


m 

8/8 

83 


CAU 

CO 

124 



30 


IKX 

8 

84 


Itt 

8 

125 E3UU 

38 

ESTAT 

31 


MTV 

C/8 

85 


8VI 

A/0EE8 

124 ESTAT: 

SS 

2 

32 


UI 

8/EXT 

84 


OP 

8 

127 



33 


MW 

8/A 

87 


22 

B0K 

128 GJU: 



34 


txx 

H 

es 


2W 

UI? 

12? CAFT: 

8S 

C 

35 


MW 

8/3 

87 




130 

88 

CSTAT 

34 


m 

H 

79 




131 CSTAT: 

BS 

•> 

37 


m 

8/C 

91 


OSES 


132 



38 




72 

SUE: 



133 90 REA0C 


37 

ffi'EK: 



73 


08 

00H/OAH/ 

'MiSTHTAX ERRO*#«'/0S8f)AH 

40 


IWI 

C/0PE8 








41 


UI 

3.08U 








42 


CAL 

ISIS 








43 


LW 

OSTAT 








44 


OJA 

A 








45 


jJG 

EXCEPT 








44 


mu 

DAFT 








47 


SHU) 

RAFT 








48 


SHLB 

CAFT 

74 


38 

'TYPE HELP FOR C0KKAKJ USTI86' 

i03H/OA8/OAH/OfE3i 

47 : 

?£>BF : 










50 


XVI 

C/PEAO 








51 


LXI 

B/R8U 








52 


CAU 

ISIS 








53 


IDA 

RSTAI 








54 


ORA 

A 
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OR 16 SWAL PA€iE IS 

QE POOR QUALITY 
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ORIGINAL PA^E jg 

OF POO/? OYAUTy 






Iw 


— J 







H. C 

. 







S 


>— » 







J r* 

z 







*JL 


tj— 







z 








H 


- 







Li O 

A 







W 

- 

-X 



h* 




x • » 

’r 






;r. 

•Itl 

UJ 

1 



L. 




H 

tL 






LI 

0L 

X 

ui 







i — 

A 

•ji 





« 


ri 

XT 



z 




c CO 


UJ 





ij 


u. 




A- 




fc- J 

}— 

Hi 





z 

•jl 


s. 



x 

• 



• 

?n 

u 





Ui 

H 

i- 

2 



1 ft 

ri 

• 


Ui 

H 

X 





X 


in 




r* 

.> 

X 


3 UJ 


UJ 





t - *: 

z 

H* 

A 




r-j 

A 


Hi j— 

-*r»- 

-A 





Ul 

o 

i 

ui 



X 


r* 


L. -J 


:*> 




S3 

•X 

r-4 


r* 




uJ 

t- 


Hi 

> 

Li 




Ua 

4> 

*-* 

til 

i.rj 



x 


X 


f“ i-L. 

x 





H 


a. 

-u. 

rA 



rH 

X 

z 


*jt 


L. 




X 

2 

o 

t— 

— j 





2 


Hi • 

£2 





p- 

u 






X 

f* 

Q 


-J 

Li 

X 





X 

j— 

Z 

in 



— J 


li 



Hi 

LJ 




a: 

Cu 

i.rt 

K 

c 




3 

z 


s a 

ii. 

T* 




3 

f— I 

r-i 


ii! 



X 

jl 

Hi 


Ui CJ 

Hi 

x 




p 

•J 

J 

C3 

!— 



H 

X 



_J •• 

*-i 





li. 



U) 

«J 



U3 

o 

-u 


S !l 

U 

u 

• 




1.1 

UJ 

UI 

l-u 





CJ 

X 


z - 

X 


iTi 



i/i 

> 


*■* 

!_r_ 



o 


z 


u z 

n 

X 

u 



UJ 

x 

r° 

> 




c-j 

• 



iC. z 


WJ 

r- 



J 

x 


• « 

r 



CN 

X 

r* 


ui •• 

u. 

:- 

X 



•2 


X 

X 

o 



=. ■* 

UJ 

}- 


X ^ 

HJ 

J 

x 



X 

«r 

hh 

w 

u. 




L— 

r-* 




ri 




u. 

« ■ 



X 



*.i 

-M 

j 

iS 


Li 2 

•A 

vL 

•3 


r— 

in 

H 

• 

VM 




A- 




X .«i 

X 




CO 

wS 


'T 


2 




u. 

X 


1 - c 

Hi 

».* 

• IU 

x 


w 

X 

«# 

X 


”1 







u» 

W 

-J 







«»\ 

r- 



X 

•— « 



Ui Z 

m 

x 




Q 

in 

C 

2 

X* 



x 

ui 

o 


Li. •-> 


Hi 

•n 



z 

X 

z 

w 

ui 




J 

•■i. 


j-i jj 

— * 




u_; 

x 

U! 

X 

J— 

j 



: /» 

w 



L 


H 




T* 

*— 


u» 

U i 



U_ 

z 

«•: 

UM 


•—* — - 

— I \u. 

Li 


z 

t 


z 


M 

w 

IT' 




Ui 



•J Q 

s 

N 

H 

i 

w 

-A 

Hi 

SJL 

— * 





»n 

Hi 

•• 


X 

•X 


»». 

u 

u- 

— 


2 




in 


Z i— — 



x 

j *• 





t.l 

UJ 

X 




c 

! A 


r.n •• 


:n 

S3 

j : • i 

KA- 



1; 

M 

x 






•• 

r ■*• 

Si 


* 

• j 

L. 

LLi 

V 

- 

H> 




u) 


:r, 


- 


•H 

: .i ; 

1 

*» 

*— 

VI_ 

— ♦ 

u 





X 


: • 
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Listing D- 20 Listing of HKLP.COM File 


920 AVAILABLE: F-'ILTEK MODULE'S : 

iiaMf i- e-: Input ft ask 
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Listing D-21 Listing of HELP . FRE File 


2920 FILTER PREPARATION PROCEDURE *• 

i. LTir'iSE KPHOM HEMOflY OF 2920 WITH 20 MINUTF EXf'05UftE IJNOtCft UO l.HiHf 


ORIGINAL PAGE IS 126 

OF POOR QUALITY 



Listing D- 22 Listing of HELP. PRO File 



Tms commands copies al.l. the help eile/s to the 
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Listing D-23 Listing of HELP. COP File 



720 PROGRAMMING EXAMPLE 


ORIGINAL PAGE IS 
OF POOR QUALITY 
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Listing D- 24 Listing of HELP . EXA File 


APPENDIX E 


SUPPORT PROGRAMS 


Appendix E contains listings and sample runs of 
three support programs used in this investigation. TRCLC 
calculates the transfer function of the analog circuit of 
Figure 1-2 as a function of a normalized value of Rl. The 
program illustrates the transfer function drift for 
various values of Rl. 

MAGCK calculates the frequency response of the 
third order Butterworth transfer function H(s) and the 
frequency response of the Bilinear Transform resultant 
H(z). The output illustrates the warping properties of 
the Bilinear Transform and the theoretical frequency 
response characteristics of the digital filter. 

MAGND calculates the magnitude response of the 
numerator and denominator of H(z) separately. The 
overshoot characteristics of the recursive and ncn- 
recursive direct form algorithms are illustrated by the 
resul ts . 
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LISTING E-l : TRCLC.BAS SOURCE CODE 


100 ! TRCLC IS A PROGRAM WHICH ALLOWS THE USER TO 

101 ! ENTER A NORMALIZED VALUE OF Rl, THE PROGRAM 

102 ! THEN CALCULATES THE TRANSFER FUNCTION OF THE 

103 I RESULTING FUNCTION. THE BASIC TRANSFER FUNCTION 

104 ! IS DERIVED FROM THE CIRCUIT DFRF USES AS A PRE- 

105 ! SAMPLING FILTER. 

106 ! 

4l0 R2=R3=C1=C2=C3=1 

111 PRINT “ Enter normalized value of Rl"; 

£ 1 2 INPUT Rl 

ul3 PRINT " Enter delta frequency " ; 

t l4 INPUT D 

15 A=R1/ ( 2*PI ) ~ 3 %\B=2*R1/ ( 2*PI ) ~ 2%\C=2*R1/ ( 2*PI ) 

V116 PRINT "FREQ- (HZ ) " , “S-TRAN (DB) " ,"S-PH (DEG) ” 
fl20 FOR N%=0% TO 25% 

130 F1=N%*D 

140 Wl=2 . *PI *F1 
150 T1=1.-B*W1~2% 

151 T2=C*W1-A*W1~3% 

160 Ml=l./SQR(Tl*Ti+T2*T2) 

170 P1=-ATN (T2/T1) 

180 IF T1<0 THEN Pl=Pl-PI *SGN (T 2 ) 

190 M 1=20 . *LOG10 (Ml ) 

200 Pl=P 1*180 . /PI 
210 PRINT FI, Ml, PI 
220 NEXT N% 

221 PRINT\PRINT\PRINT\PRINT 
230 GOTO 100 
240 END 
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Enter 

normalized value of 

Rl? 1.1 

Enter 

delta frequency? .1 


FREQ- 

(HZ) S-TRAN(DB) 

S-PH(DEG) 

0 

0 

0 

. 1 

-. 190723E-1 

-12.6162 

. 2 

-. 761015E-1 

-25.3051 

. 3 

-.172346 

-38.1642 

.4 

-.315523 

-51.3264 

. 5 

-.526579 

-64.9424 

. 6 

-.84524 

-79.1223 

. 7 

-1.32885 

-93.838 

.8 

-2.03793 

-108.825 

. 9 

-3.00909 

-123.575 

1 

-4.23246 

-137.49 

1 . 1 

-5.65374 

-150.095 

1.2 

-7.19881 

-161.173 

1.3 

-8.79901 

-170.737 

1.4 

-10.4033 

-178.934 

1.5 

-11.979 

174.038 

1.6 

-13.5077 

167.988 

1.7 

-14.9801 

162.744 

1.8 

-16.3929 

158.166 

1.9 

-17.7462 

154.14 

2 

-19.0417 

150.573 

2.1 

-20.2824 

147.391 

2.2 

-21.4712 

144.536 

2.3 

-22.6115 

141.959 

2.4 

-23.7064 

139.621 

2.5 

-24.759 

137.49 


SAMPLE RUN E-l: TRCLC.BAS for +10% DRIFT 



Enter normalized value of Rl? 1.0 
Enter delta frequency? .1 


FREQ- (HZ) 

S-TRAN(DB) 

S-PH (DEG) 

0 

0 

0 

. 1 

-. 434294E-5 

-11.4785 

. 2 

27794E-3 

-23.0782 

. 3 

-. 316485E-2 

-34.9451 

.4 

-. 177 524E-1 

-47.2648 

. 5 

-. 673338E-1 

-60.2551 

.6 

-.19804 

-74.1161 

.7 

-.483054 

-88.916 

.8 

-1.01109 

-104.432 

.9 

-1.851 

-120.066 

1 

-3.0103 

-135 

1.1 

-4.42724 

-148.535 

1.2 

-6.00536 

-160.331 

1.3 

-7.65431 

-170.389 

1.4 

-9.30925 . 

-178.901 

1.5 

-10.9309 

173.884 

1.6 

-12.4986 

167.731 

1.7 

-14.0032 

162.436 

1.8 

-15.4422 

157.839 

1.9 

-16.8166 

153.812 

2 

-18.1291 

150.255 

2.1 

-19.3835 

147.09 

2.2 

-20.5835 

144.253 

2.3 

-21.7329 

141.696 

2.4 

-22.8353 

139.377 

2.5 

-23.8942 

137.265 


SAMPLE RUN E-2: TRCLC.BAS for ZERO DRIFT 



Enter normalized value of Rl? .9 


Enter 

delta 

frequency? .1 


FREQ- 

(HZ) 

S-TRAN(DB) 

S-PH (DE) 

0 


0 

0 

. 1 


. 156593E-1 

-10.3362 

. 2 


. 627646E-1 

-20.8154 

. 3 


. 140392 

-31.6079 

.4 


. 242385 

-42.9332 

. 5 


. 349601 

-55.0691 

.6 


. 418226 

-68.3236 

.7 


. 367338 

-82.9291 

.8 


. 791 653E-1 

-98.8235 

.9 


-.564731 

-115.415 

1 


-1.61368 

-131.634 

1.1 


-3.00892 

-146.419 

1.2 


-4.62435 

-159.198 

1.3 


-6.33601 

-169.928 

1.4 


-8.05727 

-178.858 

1.5 


-9.73885 

173.686 

1.6 


-11.357 

167.401 

1.7 


-12.9028 

162.045 

1.8 


-14.3749 

157.426 

1.9 


-15.7758 

153.401 

2 


-17.1096 

149.859 

2.1 


-18.3811 

146.715 

2.2 


-19.595 

143.902 

2.3 


-20.7556 

141.37 

2.4 


-21.8673 

139.076 

2.5 


-22.9337 

136.987 


SAMPLE RUN E-3: TRCLC.BAS for -10% DRIFT 
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LISTING E-2 : MAGCK.BAS SOURCE CODE 


100 l MAG (NITUDE ) C(HEC)K CALCULATES THE EXPECTED 

101 ! VALUE OF THE S DOMAIN FUNCTION AND THE ACTUAL 

102 ! Z DOMAIN RESULTS USING THE ROOT POLE PLOTS 

103 ! DESCRIBED IN ANTONIOU. 

104 ! 

109 PRINT 

110 INPUT "Enter desired corner frequency (Hz) ";F0 

111 W0=2.*PI*F0 

112 N1=F0/10 

113 F=3. 2E6/768 

114 GOSUB 300 

115 PRINT "FREQ- (HZ ) " , " S-TRAN (DB) " , " Z-TRAN (DB) " ; 

116 PRINT " S-PH (DEG) " ,"Z-PH (DEG) " 

140 G=2**-4 

150 FOR N%=0% TO 25% 

160 F 1=N%*N 1 

161 W1=2.*PI*F1 

162 A =COS (Wl/F) 

163 B =SIN (Wl/F) 

170 X=FNR (A0,A1 ,A2, A3) 

171 Y=FNI (A0 ,A1 , A2 , A3 ) 

172 W=FNR(K0,K1,K2,K3) 

173 Z=FNI (K0,K1,K2,K3) 

180 M=SQR { (X“2 %+Y ~2% ) / (W~2%+Z~2%) ) 

181 P 1=ATN (Y/X ) 

182 IF X<0. THEN P1=P1+PI*SGN (Y) 

190 P2=ATN ( Z/W) 

191 IF W<0. THEN P2=P2+PI *SGN ( Z ) 

195 P=P1-P2 

196 IF P>PI THEN P=P-2*PI 

197 IF P<-PI THEN P=P+2*PI 

200 T1=W0“3%-2.*W0*W1~2% 

201 T2=2.*W0“2%*W1-W1"3% 

210 M1=W0~3%/SQR (T1*T1+T2*T2 ) 

211 P1=-ATN (T2/T1) 

212 IF T1<0 THEN P1=P 1 -P I *SGN (T 2 ) 

213 M1=.8*M1 

214 M =G*M 

215 M 1=20* LOG 10 (Ml) 

216 M =20 *LOG 10 (M) 

217 P1=P1*180 ./PI 

218 P =P *180. /PI 

220 PRINT FI, Ml, M, PI, P 

230 NEXT N% 

240 PRINT 

241 STOP 

250 DEF* FNR (C0,C1,C2,C3) 

251 FNR=C0*A~3%+C1*A"2%+C2*A-3.*C0*A*B~2%-C1*B‘'2%+C3 
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252 

265 

266 
267 
300 
310 
320 
330 
340 
350 
360 
370 
380 
390 
400 
410 

420 

421 

422 

423 

424 

425 

426 

427 

428 

429 

430 

431 
440 
500 


FNEND 

DEF* FNI (00,01,02,03) 

FNI=-C0*B~ 3 %+3.*C0*A~2%*B+2.*Cl*A*B+C2*B 
FNEND 

L=W0/(2*F*TAN (W0/(2*F) ) ) 

R=L*F/W0 

K0=8*R**3+8 *R** 2+4 *R+1 
Kl=-24.*R**3-8.*R**2+4.*R+3. 
K2=24.*R**3-8.*R**2-4. *R+3 . 
K3=-8.*R**3+8.*R**2-4.*R+1. 

K 3=K3/K0 

K2=K2/K0 

Kl=Kl/K0 

A0 = . 8/K0*2 ** 4 

Al=.8*3/K0*2**4 

A2=A1 

A3=A0 

K 0 = 1.0 

PRINT 

PRINT "Z-TRANSFORM COEFFICIENTS:" 

PRINT 


PRINT 

“ A ( 0 ) = 

" ; AO 


PRINT 

" A ( 1 ) = 

"?A1,"K(1)= 

" ;K1 

PRINT 

" A { 2 ) = 

°;A2,"K(2)= 

" ; K2 

PRINT 

PRINT 

* A ( 3 ) = 

" ; A3 , °K (3 ) = 

" ; K3 

PRINT 

" Rsc= " 

?R 



PRINT 

RETURN 

END 



Enter desired corner frequency (Hz) ? 300 


Z -TRANSFORM COEFFICIENTS: 


A ( 0 ) = 

. 988407E-1 




A ( 1 ) = 

.296522 

K(l) = 

-2.10202 


A ( 2 ) = 

. 296522 

K(2) = 

1.56514 


A ( 3 ) = 

. 988407E-1 

K(3) = 

-.401342 


Rsc = 

2.17266 




FREQ- 

(HZ) S-TRAN(DB) 

Z-TRAN (D3) 

S-PH (DEG) 

Z-PH (DEG) 

0 

“1.9382 

-1.9382 

0 

0 

30 

--1.9382 

-1.9382 

-11.4785 

-11.2833 

60 

-1.93848 

-1.93845 

-23.0782 

-22.6934 

90 

-1.94137 

-1.94108 

-34.9451 

-34.3812 

120 

-1.95595 

-1.95447 

-47.2648 

-46.538 

150 

-2.00553 

-2.00052 

-60.2551 

-59.3927 

180 

-2.13624 

-2.12375 

-74.1161 

-73.1691 

210 

-2.42125 

-2.39753 

-88.916 

-87.9778 

240 

-2.94929 

-2.9159 

-104.432 

-103.648 

270 

-3.7892 

-3.75945 

-120.066 

-119.607 

300 

-4.9485 

-4.9485 

-135 

-135 

330 

-6.36544 

-6.42695 

-148.535 

-149.047 

360 

-7.94356 

-8.09529 

-160.331 

-161.339 

390 

-9.59251 

-9.85655 

-170.389 

-171.84 

420 

-11.2475 

-11.6401 

-178.901 

179.261 

450 

-12.8691 

-13.4031 

173.884 

171.707 

480 

-14.4368 

-15.1232 

167.731 

165.248 

510 

-15.9414 

-16.7904 

162.436 

159.673 

540 

-17.3804 

-18.4022 

157.839 

154.814 

570 

-18.7548 

-19.9599 

153.812 

150.536 

600 

-20.0673 

-21 . 4666 

150.255 

146.737 

530 

-21 . 3217 

-22.9263 

147.09 

143.336 

660 

-22.5217 

-24.3434 

144.253 

140.267 

690 

-23.6711 

-25.722 

141.696 

137.48 

720 

-24.7735 

-27.0663 

139. 377 

134.934 

750 

-25.8324 

-28.3801 

137.265 

132.595 


SAMPLE RUN E-4: MAGCK.BAS for a 300 Hz FILTER 



Enter desired corner frequency (Hz) ? 80 


Z -TRANSFORM COEFFICIENTS: 


A ( 0 ) = 

- 249874E-2 




A ( 1 ) = 

. 749622E-2 

K(l) = 

-2.75887 


A ( 2) = 

. 749622E-2 

K(2) = 

2.54594 . 


A ( 3 ) = 

. 249874E-2 

K (3) = 

-.785511 


Rsc = 

8.27926 




FREQ- 

(HZ) S-TRAN(DB) 

Z-TRAN (OB) 

S-PH (DEG) 

Z-PH (DEG) 

0 

-1.9382 

-1.9382 

0 

0 

8 

-1.9382 

-1.9382 

-11.4785 

-11.4646 

16 

-1.93848 

-1.93848 

-23.0782 

-23.051 

24 

-1.94137 

-1.94134 

-34.9451 

-34.9052 

32 

-1.95595 

-1.95584 

-47.2648 

-47.2133 

40 

-2.00553 

-2.00517 

-60.2551 

-60.1942 

48 

-2.13624 

-2.13534 

-74.1161 

-74.0493 

56 

-2.42125 

-2.41956 

-88.916 

-88.85 

64 

-2.94929 

-2.94693 

-104.432 

-104.377 

72 

-3.7892 

-3.78712 

-120.066 

-120.034 

80 

-4.9485 

-4.9485 

-135 

-135 

88 

-6.36544 

-6.3697 

-148.535 

-148.57 

96 

-7.94356 

-7.954 

-160.331 

-160.401 

104 

-9.59251 

-9.61062 

-170.389 

-170.49 

112 

-11.2475 

-11.2743 

-178.901 

-179.029 

120 

-12.8691 

-12.9056 

173.884 

173.733 

128 

-14.4368 

-14.4835. 

167.731 

167.558 

136 

-15.9414 

-15.999 

162.436 

162.244 

144 

-17.3804 

-17.4494 

157.839 

157.629 

152 

-18.7548 

-18.8359 

153.812 

153.584 

160 

-20.0673 

-20.1611 

150.255 

150.011 

168 

-21.3217 

-21.4287 

147.09 

146.829 

176 

-22.5217 

-22.6426 

144.253 

143.977 

184 

-23.6711 

-23.8065 

141.696 

141.404 

192 

-24.7735 

-24.9241 

139. 377 

139.07 

200 

-25.8324 

-25.9986 

137.265 

136.943 


SAMPLE RUN S-5: MAGCK.BAS for a 80 Hz FILTER 



138 


Enter desired corner frequency (Hz) ? 70 


Z -TRANSFORM COEFFICIENTS: 


A ( 0 ) = 

. 169798E-2 




A ( 1 ) = 

. 509395E-2 

K(l) = 

.-2.78898 


A ( 2 ) = 

. 509395E-2 

K(2) = 

2.59964 


A ( 3 ) = 

. 169798E-2 

K(3) = 

-.809601 


Rsc = 

9.46471 




FREQ- 

(HZ) S-TRAN(DB) 

Z-TRAN (D3 ) 

S-PH (DEG) 

Z-PH (DEG) 

0 

-1.9382 

-1.9382 

0 

0 

7 

-1.9382 

-1.9382 

-11.4785 

-11.4679 

14 

-1.93848 

-1.93848 

-23.0782 

-23.0573 

21 

-1.94137 

-1.94135 

-34.9451 

-34.9145 

28 

-1.95595 

-1.95587 

-47.2648 

-47.2254 

35 

-2.00553 

-2.00526 

-60.2551 

-60.2084 

42 

-2.13624 

-2.13555 

-74.1161 

-74.065 

49 

-2.42125 

-2.41996 

-88.916 

-88.8655 

56 

-2.94929 

-2.94748 

-104.432 

-104. 39 

63 

-3.7892 

-3.78761 

-120.066 

-120 . 042 

70 

-4.9485 

-4.9485 

-135 

-135 

77 

-6.36544 

-6.3687 

-148.535 . 

-148.562 

84 

-7.94356 

-7.95155 

-160.331 

-160.384 

91 

-9.59251 

-9.60637 

-170.389 

-170.467 

98 

-11.2475 

-11.268 

-178.901 

-178.999 

105 

-12.8691 

-12.897 

173.884 

173.769 

112 

-14.4368 

-14.4726 

167.731 

167.598 

119 

-15.9414 

-15.9855 

162.436 

162.289 

126 

-17.3804 

-17.4332 

157.839 

157.678 

133 

-18.7548 

-18.8168 

153.812 

153.638 

140 

-20.0673 

-20.139 

150.255 

150.068 

147 

-21.3217 

-21.4036 

147.09 

146.89 

154 

-22.5217 

-22.6142 

144.253 

144.042 

161 

-23.6711 

-23.7746 

141.696 

141.473 

168 

-24.7735 

-24.8886 

139.377 

139.142 

175 

-25.8324 

-25.9595 

137.265 

137.018 


SAMPLE RUN E-6: MAGCK.BAS for a 70 Hz FILTER 


LISTING E-3: MAGND.BAS SOURCE CODE 


100 l MAG (NITUDE ) N (UMERATOR) & D (ENOMINATOR) CALCU- 

101 ! LATES THE PHASE AND MAGNITUDE RESPONSE OF THE 

102 ! RECURSIVE RESULTS USING THE ROOT POLE PLOTS 

103 ! DESCRIBED IN ANTONIOU. 

104 ! 

109 PRINT 

110 INPUT "Enter desired corner frequency (Hz) ”;F0 

111 W0=2.*PI*F0 

112 N1=F0/10 

113 F=3 . 2E6/768 

114 GOSUB 300 

115 PRINT " FREQ (Hz ) ” , " NUM (dB) " , "DEN (dB) " , " MAG (dB) " 

140 g=2**-4 

150 FOR N%=0 % TO 25% 

160 F 1=N %*N 1 

161 W1=2.*PI*F1 

162 A =COS (Wl/F) 

163 B =SIN(W1/F) 

170 X=FNR(A0,A1,A2,A3) 

171 Y=FNI (A0,A1,A2,A3) 

172 W=FNR(K0,K1,K2,K3) 

173 Z =FNI (K0,K1,K2,K3) 

180 M 1 =S QR ( X ~ 2 %+ Y ~ 2 % ) 

181 P1=ATN (Y/X) 

182 IF X<0. THEN P1=P1+PI*SGN (Y) 

185 M2=1/SQR(Z~2%+W~2%) 

190 P 2=ATN ( Z/W) 

191 IF W<0. THEN P2=P2+PI*SGN (Z) 

213 M1=20*LOG10 (G*M1) 

214 M2=20*LOG10 (M2) 

217 P1=P1*180./PI 

218 P2=P2*180./PI 

220 PRINT FI, Ml, M2, Ml +M 2 

230 NEXT N% 

240 PRINT 

241 STOP 

250 DEF* FNR (C0,C1,C2,C3) 

251 FNR=C0*A''3 %+Cl*A* 2%+C2*A-3 ,*C0 llt A*B"2%-Cl*B''2%+C3 

252 FNEND 

265 DEF* FNI (CO ,C1 , C 2 ,C3 ) 

266 FNI=-C0*B~3%+3.*C0*A~2%*B+2.*C1*A*B+C2*B 

267 FNEND 

300 L=W0/(2*F*TAN(W0/(2*F) ) ) 

310 R=L*F/W0 

320 K0=8*R**3+8*R**2+4*R+1 

330 Kl=-24.*R**3-8.*R**2+4.*R+3. 

340 K2=24.*R**3-8.*R**2-4.*R+3. 

350 K3=~8 . *R**3+8 . *R**2-4 . *R+1 . 


360 

K3=K3/K0 



370 

K2=K2/K0 



380 

K 1=K1/K0 



390 

A0=.8/K0*2**4 



400 

Al=.8*3/K0*2**4 



410 

A2=A 1 



420 

A3=A0 



421 

K0 = 1 . 0 



422 

PRINT , 



423 

PRINT "Z-TRANSFORM 

COEFFICIEN 

424 

PRINT 



425 

PRINT " A ( 0 ) = " ;A0 



426 

PRINT M A ( 1 ) = " ;A1, 

" K ( 1 } = 

" ; K1 

427 

PRINT " A { 2 ) = " ; A2 , 

"K (2)= 

" ; K2 

428 

PRINT " A ( 3 ) = n ; A3 , 

" K ( 3 ) = 

";K3 

429 

PRINT 



430 

PRINT M Rsc= " ; R 



431 

PRINT 



440 

RETURN 



500 

END 





Enter desired corner frequency (Hz) ? 300 


Z-TRANSFORM COEFFICIENTS: 



A ( 0 ) = . 

988407E-1 



A(l) = . 

296522 

K ( 1 ) = -2 

. 10202 

A(2)= . 

296522 

K(2)= 1 

. 56514 

A ( 3 ) = . 

988407E-1 

K(3)= -. 

401342 

Rsc= 2. 

17266 



FREQ (Hz ) 

NUM (dB) 

DEN (dB) 

MAG (dB) 

0 

-26.1219 

24.1837 

-1.9382 

30 

-26.1285 

24.1903 

-1.9382 

60 

-26.1486 

24.2101 

-1.93845 

90 

-26.1819 

24.2408 

-1.94108 

120 

-26.2287 

24.2742 

-1.95447 

150 

-26.2889 

24.2884 

-2.00052 

180 

-26.3626 

24.2388 

-2.12375 

210 

-26.4499 

24.0524 

-2.39753 

240 

-26.5509 

23.635 

-2.9159 

270 

-26.6656 

22.9062 

-3.75945 

300 

-26.7943 

21.8458 

-4.9485 

330 

-26.9369 

20. 51 

-6.42695 

360 

-27.0938 

18.9985 

-8.09529 

390 

-27.2651 

17.4085 

-9.85655 

420 

-27.4509 

15.8108 

-11.6401 

450 

-27.6514 

14.2483 

-13.4031 

480 

-27.867 

12.7439 

-15.1232 

510 

-28.0978 

11.3074 

-16.7904 

540 

-28.3441 

9.94192 

-18.4022 

570 

-28.6063 

8.64637 

-19.9599 

600 

-28.8846 

7.41794 

-21.4666 

630 

-29.1793 

6.25297 

-22.9263 

660 

-29.4909 

5.14756 

-24.3434 

690 

-29.8198 

4.09783 

-25.722 

720 

-30.1664 

3.10014 

-27.0663 

750 

-30.5312 

2.15108 

-28.3801 


SAMPLE RUN E-7: MAGND.BAS for a 300 Hz FILTER 



Enter desired corner frequency (Hz) ? 80 


Z -TRANSFORM COEFFICIENTS: 


A(0)« . 

249874E-2 



A(l)= . 

749622E-2 

K(l)= -2 

. 75887 

A ( 2 ) = . 

749622E-2 

K(2)= 2 

. 54594 

A(3)= . 

249874E-2 

K(3)= -. 

785511 

Rsc= 8. 

27926 



FREQ ( Hz ) 

NUM (dB) 

DEN (dB) 

MAG (dB) 

0 

-58.0662 

56.128 

-1.9382 

8 

-58.0667 

56.1285 

-1.9382 

16 

-58.0681 

56.1296 

-1.93848 

24 

-58.0705 

56.1291 

-1.94134 

32 

-58.0738 

56.1179 

-1.95584 

40 

-58.078 

56.0729 

-2.00517 

48 

-58.0833 

55.9479 

-2.13534 

56 

-58.0894 

55.6699 

-2.41956 

64 

-58.0965 

55.1496 

-2.94693 

72 

-58.1046 

54.3175 

-3.78712 

80 

-58.1136 

53.1651 

-4.9485 

88 

-58.1236 

51.7539 

-6.3697 

96 

-58.1345 

50.1805 

-7.954 

104 

-58.1464 

48.5358 

-9.61062 

112 

-58.1592 

46.8849 

-11.2743 

120 

-58 . 173 

45.2674 

-12.9056 

128 

-58.1877 

43.7042 

-14.4835 

136 

-58.2034 

42.2044 

-15.999 

144 

-58.2201 

40.7706 

-17.4494 

152 

-58.2377 

39.4018 

-18.8359 

160 

-58.2563 

38.0952 

-20.1611 

168 

-58. 2758 

36.8471 

-21.4287 

176 

-58.2963 

35.6537 

-22.6426 

184 

-58.3178 

34.5112 

-23.8065 

192 

-58.3402 

33.4161 

-24.9241 

200 

-58.3636 

32.365 

-25.9986 


SAMPLE RUN E-8: MAGND.BAS for a 80 Hz FILTER 



Enter desired corner frequency (Hz) ? 70 


Z -TRANSFORM COEFFICIENTS: 


A (0 ) = 

. 169798E-2 



A (1) = 

. 509395E-2 

K ( 1 ) = -2 

.78898 

A ( 2 ) =. 

. 509395E-2 

K ( 2 ) = 2 

.59964 

A ( 3 ) = 

. 169798E-2 

K(3)= -. 

809601 

Rsc = 

9.46471 



FREQ (H 

z) NUM (dB) 

DEN (dB) 

MAG (dB) 

0 

-61.4219 

59.4837 

-1.9382 

7 

-61.4223 

59.4841 

-1.9382 

14 

-61.4234 

59.4849 

-1.93848 

21 

-61.4252 

59.4838 

-1.94135 

28 

-61.4277 

59.4719 

-1.95587 

35 

-61.431 

59.4257 

-2.00526 

42 

-61.435 

59.2994 

-2.13555 

49 

' -61.4397 

59.0198 

-2.41996 

56 

-61.4452 

58.4977 

-2.94748 

63 

-61.4513 

57.6637 

-3.78761 

70 

-61.4582 

56.5097 

-4.9485 

77 

-61.4659 

55.0972 

-6.3687 

84 

-61.4742 

53.5227 

-7.95155 

91 

-61 . 4833 

51.8769 

-9.60637 

98 

-61.4931 

50.2251 

-11.268 

105 

-61.5037 

48.6067 

-12.897 

112 

-61.5149 

47.0424 

-14.4726 

119 

-61 . 527 

45.5415 

-15.9855 

126 

-61.5397 

44.1065 

-17.4332 

133 

-61.5532 

42.7364 

-18.8168 

140 

-61.5674 

41.4283 

-20.139 

147 

-61.5823 

40.1788 

-21.4036 

154 

-61.598 

38.9838 

-22.6142 

161 

-61.6144 

37.8397 

-23.7746 

168 

-61.6315 

36.7429 

-24.8886 

175 

-61.6494 

35.69 

-25.9595 


SAMPLE RUN E-9: 

MAGND.BAS for 

a 70 Hz FILTER 



GLOSSARY 


A/D 

ADC 

ALU 

ASP 

C i 

D/A 

DAC 

DAR 

DFRF 

S rm 

E rms 

EPROM 

F c 

p 

C cry 

FIR 

F pcr 

F s 

F 

c smax 
F 

r smin 
H ( s) 

H ( z) 


gain of section i 

Analog to Digital 

Analog to Digital Convertors 

Arithmetic Logic Unit 

Analog Signal Processor 

gain scale factor 

Digital to Analog 

Digital to Analog Converters 

Digital Analog Register 

Dryden Flight Research Facility 

finite multiplication accuracy error 

finite signal accuracy error 

Electrically Programmable Read Only Memory 

filter corner frequency (Hz) 

processor crystal frequency (Hz) 

Finite Impulse Response 

program cycle rate (Hz) 

filter sample frequency (Hz) 

maximum 2920 sample rate (Hz) 

minimum acceptable 2920 sample frequency (Hz) 

analog transfer function 

discrete transfer function 


144 




1 


145 


PSFSS 


X(nT) 


Y (nT) 


Input-Output 

Infinite Impulse Response filter structure 

pre-warping factor 

Number of multiplications 

Number of program steps in filter code 

Pulse Coded Modulation 

Pre-Sampling Filter 

Pre-Sampling Filter Support Software 
binary Resolution of input section i 
binary Resolution of signal 
Random Access Memory 

sample rate to corner frequency ratio (Hz/Rad) 

Laplacean operator 

binary Scale Factor for section i 

Signal Conditioning Buffers 

analog frequency variable (radians/sec) 

digital frequency variable (radians/sec) 

filter corner frequency (radians/sec) 

discrete section input at time nT 

discrete section output at time nT 

discrete time operator 

Intel 2920 Analog Signal Processor 


2920 


